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Abstract. Predictions of the concordance cosmological model (CCM) of the structures in the environment of large spiral 
galaxies are compared with observed properties of Local Group galaxies. Five new most probably irreconcilable problems 
are uncovered: 1) A great variety of published CCM models consistently predict some form of relation between dark-matter- 
mass and luminosity for the Milky Way (MW) satellite galaxies, but none is observed. 2) The mass function of luminous 
sub-halos predicted by the CCM contains too few satellites with dark matter (DM) mass a 10 7 M e within their innermost 
300 pc than in the case of the MW satellites. 3) The Local Group galaxies and data from extragalactic surveys indicate there 
is a correlation between bulge-mass and the number of luminous satellites that is not predicted by the CCM. 4) The 13 new 
ultra-faint MW satellites define a disc-of-satellites (DoS) that is virtually identical to the DoS previously found for the 11 
classical MW satellites, implying that most of the 24 MW satellites are correlated in phase-space. 5) The occurrence of two 
MW-type DM halo masses hosting MW-like galaxies is unlikely in the CCM. However, the properties of the Local Group 
galaxies provide information leading to a solution of the above problems. The DoS and bulge-satellite correlation suggest that 
dissipational events forming bulges are related to the processes forming phase-space correlated satellite populations. These 
events are well known to occur since in galaxy encounters energy and angular momentum are expelled in the form of tidal tails, 
which can fragment to form populations of tidal-dwarf galaxies (TDGs) and associated star clusters. If Local Group satellite 
galaxies are to be interpreted as TDGs then the substructure predictions of the CCM are internally in conflict. All findings 
thus suggest that the CCM does not account for the Local Group observations and that therefore existing as well as new viable 
alternatives have to be further explored. These are discussed and natural solutions for the above problems emerge. 
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1. Introduction 

Our understanding of the cosmological world relies on two fun- 
damental assumptions: 1) The validity of General Relativity, 
and 2) conservation of matter since the Big Bang. Both as- 
sumptions yield the concordance cosmological model (CCM), 
according to which an initial inflationary period is followed by 
(exotic, i.e., non-baryonic) dark-matter (DM) structures form- 
ing and then accreting baryonic matter, which fuels star for- 
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mation in the emerging galaxies, and according to which dark 
energy (represented by a cosmological constant A) drives the 
acceleration of the Universe at a later epoch. One important 
way to test assumption (1) is to compare the phase-space prop- 
erties of the nearest galaxies with the expectations of the CCM. 
These tests are the focus of the present contribution. 

The possibility of th e existence of DM was considered 



more than 85 years ago (lEinsteinl 1 1 92 It lOort 111932c IZwickv 
1933b . and has been under heav y theoretical and experimen- 



tal scrutiny (IBertone et al.l 120051) si nce the discovery of non- 
Keplerian galactic rotation curves bv lRubin&Fbrdl ( ll970l) and 
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their verification and full establishment bv lBosmal (119811) . The 
existence of DM is popularly assumed because it complies with 
the General Theory of Relativity, and therefore Newtonian dy- 
namics, in the weak-field limit. Newtonian dynamics is the 
simplest form of gravitational dynamics given that the equa- 
tions of motion are linear in the potential, and is thus readily 
accessible to numerical simulations of cosmic evolution, upon 
which the concordance scenario of structure formation is based 



dBlumenfhal et al. 1119841) . 

The concordance bottom-up scenario of structure formation 
involving the repeated accretion of clumps of cold dark matter 
(CDM) is assumed to operate throughout the Universe on all 
scales. CDM particles with masses of order of about 100 GeV 
are the preferred candidates to account for constraints placed 
on the matter density, SIm, of thermal relics with realistic cross- 



sections (see, e.g., eq. 28 of iBertone et alj|2005l) . For lighter 
particle candidates, the damping scale becomes too large: for 
instance, a hot DM (HDM) particle candidate (;7Jhdm ~ few 
eV) would have a free-streaming length of « lOOMpc lead- 
ing to too little power at the small-scale end of the matter 
power spectrum. The existence of galaxies at redshift z ~ 6 
implies that the coherence scale should have been smaller than 
100 kpc or so, meaning that warm DM (WDM) parti cles with 
mass otwdm ~ 1 - 10 keV are close to being ruled out ( Peacockl 
20031) . 

CDM is a concept that, together with the cosmological 
constant (A), has been motivated primarily by large-scale ob- 
servations of, e.g., the cosmic microwave background (CMB ) 
radiation (WMAP, ISpergel et al.l 120071: iKomatsu et alj|2009b. 
the accelerating universe jRiess et alJ ll998t iPerlmutter et al 



19991) . or the power spectru m of density perturbations from the 
SPSS dTegmark et al.l2004 and the 2dF galaxy redshift survey 
( Cole et al.ll2.005h . all of which serve as empirical benchmarks 
for calibrating and constraining theoretical scenarios and cos- 
mological models. This concordance ACDM mod el is consis- 



tent w ith observations on the Gpc to Mpc scales (Re yes et al 



201 Oh . but it implies that the Universe evolves towards an in- 
finite energy contend due to the creation of vacuu m energy 
from dark-energy-driven accelerated expansion (e.g. Peacockl 
1999pl . Less problematically perhaps, but nevertheless note- 
worthy, the DM particle cannot be contained in the Standard 
Model of particle physics witho ut implying a significant re- 
vision of particle physics (e.g. Peacockl 19991) . Strong evi- 
dence for the existence of DM has been suggested by the 
observations of the i nteracting galaxy-c luster pair 1E0657-56 



(the "Bullet cluster", IClowe et alj 120061) . The velocity of the 



sub-cluster relative to the large cluster has since been calcu- 
lated to be ab out 3000 kms 1 so that the o bserved morphol- 
ogy can arise ( Mastropietro & Burkertll2008 ). But according to 



1 One may refer to this issue as the "cosmological energy catas- 
trophy" in allusion to the black body UV catastrophy, which led 
Max Planck to heuristically introduce an auxiliary (= HilfsgroBe in 
German) number h, to reproduce the black body spectrum. 

2 Energy conservation is a problematical issue in General Relativity 
(GR). The stress-momentum-energy tensor is a pseudo tensor and so 
is not invariant under a transformation to a different coordinate system 
and back. This may perhaps be considered to indicate that GR may not 
be complete. 



Angus & McGaughl d2008l) and iLee & Komatsul d2010l) . such 
high velocities between a sub-cluster and a main galaxy cluster 
are virtually excluded in the CCM. Near the centre of lens- 
galaxies, the observed delay times between the multiple im- 
ages of strongly lensed background sources cannot be under- 
stood if the galaxy has a standard (NFW or isothermal) DM 
content and if, at the same time, the Hubble constant has a 
classical value of 70 km s _1 Mpc -1 : the solution is either to de- 
crease the Hubble constant (in disagreement with other obser- 
vations), or to consider the known baryonic matter (with con- 
stant rnassjojjigjuratio^^heone and only source of the lens- 
ing (IKochanek & Schechterl 120041) . On Local Volume scales 
(within about 8 Mpc), it has been pointed out that the Local 
Void contains far fewer dwarf galaxies than expected if the 
CCM were true. At the same time, there are too many large 
galaxies in the less crowded parts such that the arrangement of 
massive galaxies i n the Local Volume is le ss than 1 per cent 
likely in the CCM dPeebles & NusseifcOlOl) . 

This discussion highlights that there are important unsolved 
issues in the CCM. This clearly means that substantial effort 
is required to understand the problems, to perhaps distill ad- 
ditional clues from the data that can provide solutions, and to 
improve the theory. 

Galaxy formation and evolution is a process that happens 
on scales much smaller than 1 Mpc. Ironically, a major limita- 
tion of our ability to develop a physically consistent model of 
how galaxies evolved out of the dark comes from incomplete 
knowledge of the Local Group, in particular from the lack of 
understanding of the structure and distribution of dwarf satel- 
lite galaxies. But, over the past few years, a steady flow of new 
results from nearby galaxies including the Milky Way (MW) 
and the improving numerical resolution of computational stud- 
ies of galaxy formation have allowed ever more rigorous tests 
of the CCM. 

According to the DM hypothesis, galaxies must have as- 
sembled by means of accretion and numerous mergers of 
smaller DM halos. Therefore, galaxies such as the MW 
should be swarmed by hundreds to thous ands of these halos 
( Moore et al. 1999at Diemand et al. 2008 ). whereby the num- 
ber of sub-halos i s smaller in WDM than in CDM models 
dKnebe et al J2008h . Furthermore, the tnaxial nature of the flow 
of matter at formation would make it impossible to dest roy halo 
substructure by violent relaxation dBoilv et al.l 120044) . These 
sub-halos should be distributed approximately isotropically 
about their host, and have a mass function such that the number 
of sub-halos in the mass inte rval M v j r , M VIJ + dM w \ r is approxi- 



mately dN oc ATj 9 dM VlI dGao et al.ll2004l) . 

In contrast to this expectation, only a few dozen shin- 
ing satellites have been found around both the MW and 
Andromeda (M31), while the next largest disc galaxy in the 
Local Group, M33, has no known satellites. The MW hosts 
the 11 "classical" (brightest) satellites, while 13 additional 
"new" and mostly ultra-faint satellite galaxies have been dis- 
covered in the past 15 years primarily by the Sloan Digital 
Sky Survey (SDSSfl While the MW satellites are distributed 



3 For convenience, the 11 brightest satellite galaxies are here re- 
ferred to as the "classical" satellites because these were known before 
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highly anisotropically (e.g. Klimentowski et al. 20 id) , obser 



vations of the internal kinematics (velocity dispersion) of the 
satellites suggest they are the most DM do minated galaxies 



known (e.g. fig. 15 in ISimon & Gehal l2007). That is, the ve 



locity dispersions of their stars seem to be defined by an un- 
seen mass component: the stars are moving faster than can be 
accounted for by their luminous matter. The known satellites 
may therefore be the luminous "tip of the iceberg" of the vast 
number of dark sub-halos orbiting major galaxies such as the 
MW. 

Much theoretical effort has been invested in solving 
the problem that the number of luminous satellites is so 
much smaller than the number of DM-halos predicted by 
the currently favoured concordance ACDM hypothesis: stel- 
lar feedback and heating processes limit baryonic growth, 
re-ionisation stops low-mass DM halos from accreting suffi- 
cient gas to form stars, and tidal forces from the host halo 



of DM sub-halos (Dekel & Silkl Il986t bekel & Wool 120031: 


Maccio et al. 2009: KoDosovetal. 


20091 


Okamoto & Frenk 


20091 iKirbv et alj|2009l Shava et al. 


20091: 


Busha et al.ll2010t 


Maccio et al.l |20 10). This impressive and important theoreti- 



cal effort has led to a detailed quantification of the DM-mass- 
luminosity relation of MW satellite galaxies. Moreover, the dis- 
covery of new (ultra-faint) dSph satellites around the MW sug- 
gests the validity of the "tip of the iceberg" notion. These lines 
of reasoning have generally led to the understanding that within 
the AC DM cosmology, no serious small-sca le issues are appar- 
ent fe.g. lTollerud et al.ll2008rJPrimackll2009b . 



In this contribution we test whether the CCM can be viewed 
as a correct description of the Universe by studying generic 
properties of the Local Grourfl which is a typical environ- 
ment for galaxies - the Local Group properties must conform 
to the CCM if it is to be valid universally. To test this hy- 
pothesis, we critically examine state-of-the art models calcu- 
lated within the CDM and WDM framework by a number of 
independent research groups developed to explain the proper- 
ties of the faint satellite galaxies, by comparing them with the 
following observations: the mass-luminosity relation for dSph 
satellites of the Milky Way (Sect. the mass-distribution of 
luminous-satellite halo-masses (Sect. [3}; and the observed re- 
lation between the bulge mass of the host galaxy and the num- 
ber of satellites (Sect. |4j. The question of whether the Disc- 
of-Satellites (DoS) exists, and if in fact the latest MW satellite 
discoveries follow the DoS, or whether the existence of the DoS 
is challenged by them, is addressed in Sect. [5] In Sect. |5j the 
observed invariance of late-type baryonic galaxies is also dis- 
cussed in the context of the Local Group. In these sections it 
emerges that the CCM has problems relating to the observed 
data. In Sect.|6]the problems are interpreted as clues to a possi- 



the SDSS era. These include the LMC and the SMC with the others 
being dwarf spheroidals. The other, more recently discovered satellites 
are fainter than the faintest "classical" satellites (UMi and Draco), and 
these are called the "new" or the "ultra-faint" satellites or dwarfs (see 
Table©. 

Useful reviews of the Local Group are provided bv lMated l fl998l) 
and Ivan den Berghl dl999h . 



ble solution of the origin of the satellite galaxies. The implica- 
tions of testing the CCM on the Local Group for gravitational 
theories are also discussed. Conclusions regarding the conse- 
quences of this are drawn in Sect. Q 

2. The satellite mass - luminosity relation 
(problem 1) 

Our understanding of the physical world relies on some fun- 
damental physical principles. Among them is the conservation 
of energy. This concept implies that it is increasingly more dif- 
ficult to unbind sub-components from a host system with in- 
creasing host binding energy. 

Within the DM hypothesis, the principle of energy conser- 
vation therefore governs how DM potentials fill-up with matter. 
There are two broadly different physical models exploring the 
consequences of this, namely models of DM halos based on 
internal energy sources (mostly stellar feedback), and models 
based on external energy input (mostly ionisation radiation). In 
the following, the observational mass-luminosity data for the 
known satellite galaxies are discussed, and the data are then 
compared to the theoretical results that are calculated within 
the CCM. 

2.1. The observational data 

Based on high quality measurements of i ndividual stellar line- 
of-sight velocities in the satellite galaxies. IStrieari et alJ(l2008l) 
(hereinafter S08) calculate dynamical masses, Mo.3kp C , within 
the inner 0.3 kpc of 18 MW dSph satellite galaxies over a wide 
range of luminosities (10 3 < L/L Q < 10 7 ). The LMC and SMC 
are excluded, as is Sagittarius because it is currently experienc- 
ing significant tidal disturbance. S08 significantly improve the 
previous works by using larger stellar data sets and more than 
double the number of dwarf galaxies, and by applying more 
detailed m ass modelling . Their r esults confirm t he earlier sug 



gestion by [Mateo et al. I d 1993b. iMateol d 1998b . iGilmore et al 



(2007|), and iPenarrubia et al.l (120081) that the satellites share a 
common DM mass scale of about 10 7 M , "and conclusively 
establish" (S08) this common mass scale. 

The finding of S08 can be quantified by writing 



log I0 M .3 kpc 

= log 10 M + Klog 10 L, 



(1) 



and by evaluating the slope, k, and the scaling, Mq. S08 derive 
k = 0.03 + 0. 03 andMp « 10 7 M o . Using the Dexter Java ap- 



plication of iDemleitner et al.l (120011) . a nonlinear, asymmetric 
error weighted least squares fit to the S08 measurements repro- 
duces the common mass and slope found by S08, as can be seen 
from the parameters listed in TableQ] By excluding the least lu- 
minous dSph data point, one obtains the same result (Table[TJ. 
It follows from Eq. Qjthat 

(M . 3kpc ) I/K 



M X J K L 



(k * 0), 

(K = 0). 



M .3kpc = M (* = 0). (2) 

This central mass of the DM halo can be tied by means of high- 
resolution CDM simulations to the total halo virial mass before 
its fall into the host halo (S08, see also Sect. [3), 

M vir = (Mo. 3k pc) 1/0 - 35 X l0~ n M Q , (3) 
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yielding M vil = 10 9 M G for M .3kpc = 10 7 M o (the common- 
mass scale for k = 0). Thus, substituting Mn.3kp C into Eq. [3] 
using Eq. |2 with k + 0, leads to 



(M vll -) a35/K = M ( ! /K xl0- (llx(O5 ^L. 



(4) 



This value of the halo mass near 10 M B for the satellites in the 



S08 sa mple is confirmed by a new analysis, in which lWolf et al 
(2010) show that the mass can be derived from a velocity dis- 
persion profile within the deprojected 3D half light profile with 
minimal assumptions about the velocity anisotropy. In this way 
they obtain a robust mass estimator. 

The observed 5<x lower value for 0.35 Ik = rj is thus r\ = 
2.06 (with k = 0.02 + 5 x 0.03 from TableQ}. 

2.2. Model type A: Internal energy sources 

Dekel & Silkl dl986l) and bekel & Wool d2003l) studied models 
according to which star formation in DM halos below a to- 
tal halo mass of M v ; r « 1O 12 M is governed by the thermal 
properties of the inflowing gas, which is regulated primarily by 
supernova feedback. These models demonstrate that the mass- 
to-light ratio of sub -halos follows M wu / L oc Zr 2 ^ 5 (eq. 24 o f 
Dekel & Wool 120031: see also eq. 33 of bekel & SUkTll986h . 
This app roximately fi ts the observed trend for dSph satellite 
galaxies dMateoll 19981) . 

These models thus imply that 



(A^vir) = £ L, 



(5) 



where L is the total luminosity, M v j r is the virial DM halo mass, 
77th = 5/3, and £ is a proportionality factor. In essence, this 
relation states that more-massive halos have a larger binding 
energy such that it becomes more difficult to remove matter 
from them than from less massive halos. 

Comparing with Eq. |4] and with its resulting rj value as 
given at the end of Sect. 12.11 it follows that the observed 5<x 
lower value for rj = 0.35 /k = 2.06 is in conflict with Eq. [5] 
where tfa = 5/3 = 1.67. 

2.3. Model type B1, B2: External energy source 

Bush a et al ] d2010h follow a different line of argument to ex- 
plain the dSph satellite population by employing the DM halo 
distribution from the via Lactea simulation. Here the notion is 
that re-ionisation would have affected DM halos variably, be- 
cause of an inhomogeneous matter distribution. A given DM 
halo must grow above a critical mass before re-ionisation to 
form stars or ac crete baryons. Thu s the inhomogeneous re- 
ionisation model ( Busha et al. 2010i their fig. 6) implies, upon 
extraction of the theoretical data and using the same fitting 
method as above, theoretical /(-values of 0.15-0.17. These dis- 
agree however, with the observational value of 0.02 with a sig- 
nificance of more than 4 <x, i.e. the hypothesis that the obser- 
vational data are consistent with the models can be discarded 
wit h a confidence of 9 9.99 per cent (Table [J}. 

Busha et al. ( 2010h suggest that adding scatter into the the- 



oretical description of how DM halos are filled with luminous 
baryons would reduce the discrepancy, but it is difficult to see 



Table 1. The slope of the DM-mass-luminosity relation of 
dSph satellite galaxies. Fitted parameters for Eq.Q] 



data 


K 




radius 
[pc] 


Mo 
[10 7 M o ] 


Observational: 










1 


+0.02 h 


:0.03 


300 


1.02 ±0.39 


2 


+0.02 d 


:0.03 


300 


1.01+0.40 


3 


+0.01 H 


:0.03 


300 


1.09 + 0.44 


*4 


-0.03 h 


:0.05 


600 


6.9 ± 4.9 


DM Models: 










A: feedback 


0.21 


300 




Bl: re-ionisation, SPS 


0.15 ± 


0.02 


300 


0.24 ± 0.06 


B2: re-ionisation 


0.17 ± 


0.01 


300 


0.18 + 0.02 


C: SAM 


0.42 ± 


0.02 


300 


2.0 + 0.9 


*D: Aq-D-HR 


0.17 ± 


0.02 


600 


0.41 + 0.14 


El: lkeV(WDM) 


0.23 ± 


0.04 


300 


0.069 ± 0.045 


E2: 5keV(WDM) 


0.12 ± 


0.02 


300 


0.43 ± 0.081 


F: Aq-infall 


0.13 ± 


0.01 


300 


0.32 ± 0.022 



Notes to the table: Fits to k = 0.35/?7: data 1-4 are observational val- 
ues, data A-F are models (see Sect. [2jl. Notes: 1: our fit to S08 (who 
give central 300 pc masses, 18 satellites, their fig. 1). 2: our fit to S08 
without Seg.l (faintest satellite, i.e. 17 satellites, their fig. 1). 3: our 
fit to S08 without Seg. 1 and without Hercules (i.e. 16 satellites, their 
fig. 1) . 4: our fit to the observational data plotted bv lOkamoto & Frenkl 
d2009l) ( who give cen t ral 600 pc masses, only 8 satellites, their 
fig. 1). A: lDekel & Sn3dl986h:lDekel & WodJ2003l) . stellar feedback 



(EqJSL 
to 



B 1 : our fit to Busha et al 



'2010|), 



their SPS model. B2: our fit 



etai](l201fj|). inhomogeneous re-ionisation model. C: our fit 
;t alJ l 

El: our fit to the 1 keV WDM model of IMaccio & Fontanotl d201( 



Maccio et ; 



to 



L v > 3 x 10 s L 



semi- analytical modelling (SA M), fit is for 
D: our fit to lOkamoto & Frenkl J2009I) (Aq-D-HR) . 



E2: our fit to the 5 keV WDM model o flMaccio & Fontanodd20ld). F 
our fit to the Aquarius sub-halo-infall models of lcooper* et al. 1 2010h . 
*: the entries with an asterisk are for the central 600 pc radius region. 



how this can be done without violating the actual scatter in the 
observed Mojkpc - L relation. 

2.4. Model type C: Semi-analytical modelling (SAM) 

Filling the multitude of DM halos with ba ryons given the abov e 
combined processes was investigated by IMaccio et al. I d2010l) . 
They semi-analytically modelled (SAM) DM sub-halos based 
on A^-body merger tree calculations and high-resolution re- 
computations. The authors state "We conclude that the num- 
ber and luminosity of Milky Way satellites can be naturally 
accounted for within the (A)Cold Dark Matter paradigm, and 
this should no longer be considered a problem." 

Their theoretical mass-luminosity data are plotted in their 
fig. 5, and a fit to the redshift z = data for L v > 3 x 10 s Ly.o 
satellites is listed in Table Q] The theoretical SAM data set 
shows a steep behaviour, k = 0.42. Given the observational 
data, this model is ruled out with a confidence of more than 
ten cr. 
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2.5. Model type D: High-resolution baryonic physics 
simulations (Aq-D-HR) 

The satellite population formed in a high-resolution N-body 
ACDM re-simulation with baryonic physics of one of the MW - 
type "Aquarius" halos is studied by lOkamoto & Frenk d2009l) . 
The treatment of baryonic processes include time-evolving 
photoionisation, metallicity-dependent gas cooling and photo- 
heating, supernova (SN) feedback, and chemical enrichment by 
means of SN la and II and AGB stars. Re-ionisation is included 
and the galactic winds driven by stellar feedback are assumed 
to have velocities proportional to the local velocity dispersion 
of the dark-matter halo. In these models 100 per cent of the 
SNII energy is deposited as thermal energy. Galactic winds are 
thus produced even for the least-massive dwarf galaxies. Winds 
are observed in strong starbursts induced through interactions 
rather than in self -regul ated dwarf gala xies, which may pose a 
problem for this ansatz dOtt et al.l2005l). T he details of the sim- 
ulations are provided by lOkamoto et al. (2010). The resultant 
sub-halo population with stars can, as claimed by the authors, 
reproduce the S08 common-mass scale. 

Following the same procedure as for the above models, this 
claim is tested by obtaining k from their fig. 1 (upper panel, 
red asterisks) and comparing it to the observational data also 
plotted in their fig. 1 (note that Okamoto & Frenkl 20091 plot 
the masses within 600 pc rather than 300 pc as used above). 
From their plot of the observational data, which only includes 
central-600 pc masses for the eight most luminous satellites, 
it follows that K bs,OF = -0.03 + 0.05. This is nicely consis- 
tent with the full S08 sample (18 satellites) discussed above. 
However, for their model data one finds that k = 0.17+0.02, i.e. 
the model can be discarded with a confidence of 3<x or 99.7 per 
cent. 



2.6. Model type E1, E2: WDM 



Mac cio & Fontanot d2010h present theoretical distributions of 
satellite galaxies around a MW-type host halo for different cos- 
mological models, namely ACDM and WDM with three pos- 
sible DM-particle masses of mwdm = T 2, and 5 keV. They 
perform numerical structure formation simulations and apply 
semi-analytic modelling to associate the DM sub-halos with 
luminous satellites. They suggest the luminosity function and 
mass-luminosity data of observed satellites is reproduced by 
the WDM models implying a possible lower limit to the WDM 
particle of towdm ~ 1 keV. 

The model and observational mass-luminosity data are 
compared in their fig. 5 for wwdm = 1 and 5 keV. The slopes of 
these model data are listed in Table Q] From Table [JJ it follows 
that the WDM model with towdm ~ 1 keV is ruled out with 
very high confidence (4cr or 99.99 per cent), and also has too 
few satellites fainter than My « -8 (their fig. 4). WDM models 
with otwdm ~ 5 keV are excluded at least with a 3<x or 99.7 per 
cent confidence, and, as is evident from their fig. 4, the models 
contain significantly too few satellites brighter than My = —11, 



2.7. Model type F: Infalling and disrupting dark-matter 
satellite galaxies 

Cooper et al. 1 (120101) study CDM model satellites in individual 
numerical models of dark matter halos computed within the 
Aquarius project. Semi-analytical modelling is employed to fill 
the sub-halos with visible matter, and the orbits of the infalling 
satellites are followed. General agreement with the observed 
satellites is claimed. 

Much as the other models above, in this numerical CDM 
model of substructure and satellite formation in a MW type 
host halo, the MW sub-halos fall-in stochastically and therefore 
do not agree with the observed phase-space correlated satel- 
lites, i.e. with the existence of a rotating DoS (Sect. [5] below). 
Furthermore, the presented model mass-luminosity data (their 
fig. 5) lead to a too steep slope (Table[TJ compared to the obser- 
vations and the DM-based model is excluded wit h a confidence 
of at least 99.7 per cent. In addition, fig. 5 of ICooper et al 



(1201 Oh shows a significant increase in the number of model 
satellites with a similar brightness as the faintest known satel- 
lite (Segue 1, hereinafter Seg. 1). This is in contradiction with 
the failure to find any additional satellites of this luminosity 
in the most recent data mining of the existing northern SDSS 
data, as discussed in Sect. l6.2l below. Indeed, observations sug- 
gest that Seg. 1 is a star clus ter rather than a satellite galaxy 
( Niederste-Ostholt et alj|2009l) . worsening this problem. 



2.8. Discussion 

In Fig. [TJ the latest theoretical ansatzes A-F to solve the cos- 
mological substructure problem are compared with the latest 
observational limit on the slope k of the DM-mass-luminosity 
relation of dSph satellite galaxies (Eq.Q]). 

The theoretical results always lead to a trend of luminosity 
with halo mass as a result of energy conservation. But the ob- 
served satellites do not show an increasing t re nd of luminosity 
with D M m ass, according to Mateo ( 1998b . IPenarrubia et al 



d2008l) . and IStrigari et al.l d2008l) . From Fig. Q] we note that 
seven ACDM models of the satellites deviate Act or more from 
the data, while only one (the WDM model E2 with wzwdm = 
5 keV, Table [TJ deviates more than 3cr from the data. The like- 
lihood |f| that any of the DM models describes the data is thus 
less than 0.3 per cent. 

As a caveat, the observed absence of a DM-mass- 
luminosity relation partially depends on the data for the ultra- 
faint dwarf s : inde ed, for the classical (most luminous) dSphs, 
Serra et al. ( 20091) argue that there may be a trend, k > 0, es- 
sentially because of their proposed increase in the mass of the 
Fornax dSph satellite. It is on the other hand plausible that the 
ultra-faint dwarfs do not possess any dark halo (see Sect. |6), 
and that the enclosed mass derived is due to observational ar- 
tifacts. In that case they should not be used as a possible im- 
provement for the missing satellite problem. This, however, 
would pose a problem for the DM hypothesis. 



5 The likelihood =1 -(confidence in per cent)/100 gives an indica- 
tion of how well the data can be accounted for by a given model. 
The confidence, as used throughout this text, is the probability level at 
which a model can be discarded. 
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DM model 

Fig. 1. The slope of the mass-luminosity relation, k (Eq. [D, 
for the models listed in Table Q] The observational constraints 
with confidence intervals are depicted as hatched regions (1,2, 
and 3<t region). Satellites with a larger dark-matter mass are on 
average more luminous such that the mass-luminosity relation 
has k > 0. However, the observational constraints lie in the 
region k « (see Table [TJ. The hypothesis that the data are 
consistent with any one of the models can be discarded with 
very high (at least 3cr, or more than 99.7 per cent) confidence. 



Aden et al. d2009bl) suggest that for the Hercules dSph 



satellite inter-loper stars need to be removed from the obser- 
vational sample, which would require a revision of the mass 
within 300 pc to the value Mo.3kpc 



1.9^ x 1O 6 M (instead 



of the value M . 3 kpc = 7.2™J x 10" M derived by S08). This 
new mass measurement, however, now lies more than 4 cr away 
from all ACDM-models considered above (Table [1]). Hercules 
can thus no t be und erstood in terms of a DM-dominated model. 
Ade n et al. d2009bh do state that DM-free models cannot be 
excluded (note also Fig. [6] below), or that Hercules may be 
experiencing tidal disturbances in its outer parts. Tidal distur- 
bance, however, would have to be very significant for its inner 
structure to be affected, because if one would require confor- 
mity with the theoretical DM-models its Mo.3kpc mass would 
have to have been much higher and similar to the value derived 
by S08 (~ 10 7 M Q ). Given the current Galactocentric distance 
of Hercules of 130 kpc and the result that the inner region of 
a satellite is only affected by tides after significan t tidal de- 
struction of its outer parts (iKazantzidis et al.l 120041) . this sce- 
nario is physically implausible. There are therefore three pos- 
sibilities: (i) Hercules is a DM-dominated satellite. This, how- 
ever, then implies that no logically consistent solution within 
the CDM framework is possible because its mass-luminosity 
datum would lie well away from the theoretical expectation, 
(ii) Hercules has no DM. This implies that it cannot be used 
in the mass-luminosity data analysis above and would also 
imply there to exist an additional type of DM-free satellites, 
which, however, share virtually all observable physical charac- 



teristics with the putatively DM filled satellites, (iii) Hercules 
has been significantly affected by tides. This case is physically 
implausible because of its large distance, but it would imply 
that Hercules cannot be used in the mass-luminosity analysis 
above (just as Sagittarius is excluded because of the signifi- 
cant tidal effects it is experiencing). Omitting Hercules from 
the data leads to a revised observational slope k = 0.01 + 0.03 
such that none of the conclusions reached above about the per- 
formance of the DM-models are affected. 

A point of contention for DM models of dSph satellite 
galaxies is that the DM halos grow at various rates and are also 
truncated variously due to tidal influence. The highly complex 
interplay between dark-matter accretion and orbit-induced ac- 
cretion truncation leads to the power-law mass function of DM 
halos, and at the same time would imply that the outcome in 
which all luminous DM sub-halos end up having the same DM 
mass were incompatible with the DM-theoretical expectations 
(see Sect. 01. 

Summarirising Sect. [2] while the theoretical results always 
lead to a trend of luminosity with halo mass, the observed satel- 
lites do not show this trend. The hypothesis that the CCM ac- 
counts for the data can be discarded with more than 99.7 per 
cent significance. 

3. The mass function of CDM halo masses 
(problem 2) 

One of the predictions of the ACDM hypothesis is the 
self-similarity of DM-halos down to (at least) the mass 
range of dwarf galaxies, i.e. that massive halos contain sub- 
halos of lowe r mass, w ith the same structure in a sta- 
tistical sense (Moore et ; 



for 



major review see 



all Il999al 

Del Popolo & Yesilvurd2007i The mass function of these sub- 
halos is, up to a critical mass M cr ; t , well approximated by 



ftub(M vir ) = 



dN 
dM vi , 



oc M. 



■1.9 



(6) 



where dN is the number of sub-h alos in the mass interval 
M vlI ,M vir + dM vir dGao etalj|2004 . M crit is given by M vil « 
O.OlMh with M h being the virial mass of the hosting CDM- 
halo. The virial mass, M v j, , is defined by 



4k , 
M vil = yA vir por vir , 



(7) 



where po is the critical density of the Universe and A v ;, is a 
factor such that A v i r po is the critical density at which matter 
collapses into a virialised halo, despite the overall expansion of 
the Universe. The virial radius r v j r is thereby determined by the 
density profile of the collap sed CDM-halo. F or M vil > 0.01 M h , 
the mass function steepens ( Gao et aDl2004l) . so that it is effec- 
tively cut off at a mass M max (see Eq.[8]below). It is reasonable 
to identify M max with the mass of the most massive sub-halo, 
which must be higher than M cl -i t , where the mass function be- 
gins to deviate from Eq.[6]and lower than Mh, the mass of the 
host-halo. Therefore, M cr - it < M max < Mf,. 

Thus, a halo with M vil « 10 12 M o , like the one that is 
thought to be the host of the MW, should have a population 
of sub-halos spanning several orders of magnitude in mass. 
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It is well known that, in consequence, a steep sub-halo mass 
function such as Eq. [6]predicts many more low-m ass sub-halos 



Eq.QTJand introducing the concentration parameter c = r v j,-/r s 
leads to 



than the number of observ ed faint MW satellites ( Mo ore et al 



1999atlKlvpin et al.ll 19991) . a finding commonly referred to as M(r) 



4np d c r 



the missing satellite problem. Efforts to solve this problem rely 
on physical processes that can either clear CDM-halos of all 
baryons or inhibit their gathering in the first place, which would 



+ c r 



In 1 + 



1 



The parameter 6 C can be expressed in terms of c, 



affect low-ma ss halos preferentially (e.g. [Moore et al ]| 2006l 



Li et alJ[20Iol: Sect 01. More specifically, \u et all J2010h find 



3 ln(l +c)-c/(l +c) 



(12) 



(13) 



that the mass function of luminous halos, 
essentially be flat for 1O 7 M < M vil < 



^lum(Afvir), WOUld 

10 9 M o . All sub- 



halos with My 



10 M Q would keep baryons and therefore 



^ium(A^vii-) = £sub(^vir) in this mass range. Thus, the mass func- 
tion of luminous sub-halos can be written as 



&im(A*w) = kkiM 



(8) 



with 

C*2 



0, 
1.' 



k 2 



l. 



h (io v ) 



i9^cr 2 -ai 



10 7 < if*. < 10 9 
in 9 < M™. < m 



as can be verified by setting r = r vil - in Eq.Q~2]and substituting 
M(r vir ) = M vir byEq.|7] 

If the halo is luminous, it is evident that M(r) is smaller 
than the total mass included within r, M r . However, assuming 
that the MW satellites are in virial equilibrium and that their 
dynamics is Newtonian in the weak-field limit, the mass-to- 
light ratios calculated for them are generally high and imply 
that they are DM-dominated and thus, M{r) = M r would be 
a good approximation. This relation is therefore adopted for 
the present discussion. In this approximation M(r = 0.3kpc) = 

A^0.3kpc- 



where the factors k, ensure that £ v ; r (M v i r ) is continuous where 
the power changes and A: is a normalisation constant chosen 
such that 



In principle, the paramet ers pn ( Navarro et al, 119971) . c 



JlO 7 



dBullock et alj|200lh . and A vir dMainini et al.ll2003l) depend on 
the redshift z but for the purpose of the present paper only z = 
needs to be considered, as this is valid for the local Universe. 
Thus, 



^vir(Afvir) dM vh = 1. 



(9) 



Po 



3Hl 
8ttG' 



where the Hubble 



(Soergeletal. 2007) 



(Mai nini et aDl2003h . and 



constant 

Avir - 



H 
98 



From a mathematical point of view, Eq. [8] is the probability 
distribution of luminous sub-halos. We note that the lu minous 
sub-halo mass fu nction proposed in lMoore et al. I (l2006l) is sim- 
ilar to the one in Li et al. ( 2010l) . In the high-mass part, it has 

the same slope as the mass function f or all sub-halos and flat- / ^ 

tens in the low-mass part (cf. fig. 3 in iMoore et al.ll2006l) . The l°gio© = 2.31 - 0. 109 log 10 ( | , 



lower mass limit for luminous halos i s however suggested to be 



(14) 

71kms _1 Mpc _I 
for ACDM-cosmology 



(15) 



M vir * lO s M in 
sub-halos has a\ 



Moore et al.l (120061). The ma ss function of all 



1.9 dGaoet al.ll2004 



where c is the expectation value of c as a function of M v j, . Thus, 
c decreases slowly with M v j r , while the scatter in the actual c is 
rather large, being 



3.1. NFW halos 

It is well established that the theoretical density profiles of 
galaxy-si zed CDM-halos are s imilar to a universal law, as pro- 
posed bv lNavarro et al.l (11997I) . The NFW profile is given by 



o-i ORm c = 0.174 



(16) 



(Ma ccio et al. 1120071) . The only caveat here is that the NFW 



PnfwM = 



S c po 



r/r s (1 + r/r 8 ) 



2 ■ 



(10) 



where r is the distance from the centre of the halo and po is the 
critical density of the Universe, while the characteristic radius 
r s and 6 C are mass-dependent parameters. 

By integrating pnfwW over a volume, the total mass of 
CDM within this volume is obtained. Thus, 



M(r) 



= f p(r')4nr' 2 dr' 
Jo 



(11) 



is the mass of CDM contained within a sphere of radius r 
around the centre of the CDM-halo, and M{r) = M vrr for 
r = r v ; r . Performing the integration on the right-hand side of 



profile is used t o integrate the mass, while the now-preferred 
Einasto profile ( Navarro et al. 2010i Sect. [TJ makes only a 
small difference in the central parts. 

3.2. Probing the ACDM hypothesis with M .3kpc 

S08 use the stellar motions in 18 MW satellites to calcule 
their mass within the central 300 pc, Mo.3k pc - They assume 
the satellites to be in virial equilibrium and that Newtonian 
dynamics can be applied to them. The sample from S08 can 
be enlarged to 20 satellites by including the Large Magellanic 
Cloud (LMC) and the Sma ll Magellanic Cloud (SMC), since 
van der Marel et al.l (120021) estimated the mass of the LMC 
within the innermost 8.9 kpc, Mlmc, using the same assump- 
tions as S08. This implies that M LM c = (8.7 + 4.3)xlO 9 M , of 
which the major part would have to be DM. Equations . 171 [T2l [L3l 
and Q3] have been used to create tabulated expectation values 
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Mass within 300 pc [ solar units ] 

Fig. 2. The mass function of luminous satellite problem. The 
cumulative distribution function for the mass within the central 
300 pc, Mo.3kpc> of the MW satellites (solid line) and the cumu- 
lative distribution function for Mojkpc of a sample of 10 6 CDM- 
halos picked from the parent distribution of luminous sub-halos 
(Eq.[8] dashed line). The null hypothesis is that the MW satel- 
lite M().3k P c masses are drawn from this parent distribution. The 
maximum distance between the two curves is 0.333 so that the 
null hypothesis can be discarded with 98.9 per cent confidence. 

of M(r) for NFW-halos with different M wu and it can thereby 
be seen that for a typical NFW-halo with M(r = 8.9 kpc) = 
8.7 x 1O 9 M , M(r = 0.3 kpc) = 2.13 x 1O 7 M = M . 3 kpc, and 
Mvw = 1.2xlO n M . We note that the SMC ha s about l/10thof 



the mass of the LMC ( Kallivavalil et al.ll2006h . hence the virial 
mass of its halo can be estimated as M v - n - = 1.2 x 1O 1O M , 
corresponding to Mo.3kpc = 1-51 x 1O 7 M . 

To test the shape of the MW satellite distribution function 
against the shape of the distribution of the Mo.3kpc values of 
the MW-satellites, artificial samples of 10 6 Mojkpc masses are 
generated in concordance with the ACDM hypothesis, using 
Monte Carlo simulations. As noted in Sect. 13. II Mojkpc is well 
approximated by M(r = 0.3kpc) in a CDM-dominated galaxy. 
M(r = 0.3kpc) can be calculated if M v i r and c are given, and 
the expectation value for c is a function of M vn -. The first step is 
therefore to choose a value for M v j r using uniform random de- 
viates and the probability distrib ution of luminous halos given 
in Eq. [8] (see e.g. chapter 7.2 in lPress et al.lll992l for details). 
The next step is to attribute a value for log 10 (c) to the cho- 
sen Myir. This is done by multiplying Eq. [16] with a Gaussian 
random deviate and adding the result to the value for log 10 (c), 
which is calculated from Eq.[l5] After transforming log 10 (c) to 
c, Mo.3kpc = M{r = 0.3kpc) of the given halo can be calculated 
from Eq. Q2] using Eqs. [7] and [13] These steps are repeated, 
until a sample of 10 6 Mo.3kpc values is generated. 

If two samples are given, the maximum distance between 
their cumulative distribution functions, D, can be calculated. 



Performing the KS-test, this quantity D allows an estimate of 
how likely it is that they are drawn from the same distribu- 
tion function. The null hypothesis is that the observed satel- 
lite galaxies are drawn from the theoretically calculated mass 
function of luminous halos; the parent distribution is thus as- 
sumed to be the mass function of M(0.3kpc) values of lumi- 
nous sub-halos according to the ACDM hypothesis. Assuming 
in Eq. [8] that M max = 1O H M , which is approximately the 
mass estimated for the CDM halo of the LMC, and taking 
M min = 10 7 M , leads to D = 0.333. According to the KS-test, 
given the parent distribution the probability of an even larger 
distance is 0.011. This means that the null hypothesis can be 
excluded with 98.9 per cent confidence. Both cumulative dis- 
tributions are shown in Fig. 

Omitting the LMC and SMC from the observational sam- 
ple but keeping M min = 1O 7 M and M max = 10" M in the 
theoretical sample yields D = 0.294, leading to the exclusion 
of the null hypothesis with a confidence of 95.5 per cent. In 
addition setting M max = 4 x 1O 1O M , which is the M vir that 
corresponds to the most massive Mo.3k P c in the S08 sample (i.e. 
the most massive remaining sub-halo), yields D = 0.301 lead- 
ing to exclusion of the null hypothesis with a confidence of 
96.3 per cent. The latter two tests comprise a homogeneous 
mass-sample of observed satellites as compiled by S08. 

That the mass function is expected to steepen at M cr ; t = 
0.01 Mh even increases the discrepancy between the ACDM 
hypothesis and the observations. Reinstating the LMC and 
SMC back into the observational sample and cutting off 
&ub(M vir ) at M max = 1O 1O M (with M min = 10 7 M o ), which 
would be close to M clit for the CDM-halo of the MW (see 
Sect. 0), and one order of magnitude below the estimated mass 
of the CDM-halo of the LMC, implies that D = 0.359 and an 
exclusion with 99.5 per cent confidence. 

On the other hand, setting M max = 1O 12 M (with M min = 
1O 7 M ) leads to D = 0.329 and an exclusion with 98.8 per 
cent confidence. Any reasonable uncertainty in the actual value 
of M max can therefore be excluded as an explanation of the dis- 
crepancy between the observed sample of Mojkpc and a sample 
generated based on the ACDM hypothesis. As a consequence, 
the same is true for the uncertainty in the actual mass of the 
halo of the MW, M h , since M max is linked to M h (see Sect. [3). 

Thus M max is kept at 10" M in the following. Adjusting 
the lower limit of ^i um (M viT ) from 10 7 M Q to 10 8 M Q then leads 
to D - 0.319 and an exclusion of the null-hypothesis with a 
confidence of 98.4 per cent. The mass of 10 8 M o is the M vil - 
suggested by the lowest Mo.3kpc in the sample from S08. We 
note that the likelihood decreases with decreasing M max , be- 
cause of the overabundance of Mq jk pc ~ 10 7 M halos becom- 
ing more prominent in the observational sample. 



6 Monte Carlo experiments are used to quantify the confidence val- 
ues for the KS-tests: Drawing the corresponding number of sub-halo 
masses (e.g. 20 as in this case) from Eq.[8] D' is calculated. This is 
repeated 10 5 times. Counting of D' values gives the fraction of cases 
when D' > D, where D is the actually obtained D' value from the 
data (e.g. D = 0.333 in this case). These fractions are reported here as 
likelihood values, and are about half as large as t he probability valu es 
obtained using approximate methods, as, e.g., bv lPress etal.l(l992h . 
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S08 suggest that £i U mC^vir) might even be cut off below a 
mass of as 1O 9 M , either because halos below that mass do 
not contain baryons or do not form at all. Indeed, modifying 
^ium(A^vir) given by Eq. [8] accordingly, results in an agreement 
between the theoretical distribution and the data (D = 0.188 
with an exclusion confidence of only 70 per cent). A ^i um (M v ; r ) 
with a lower mass limit of 1O 9 M is however in disagree- 
ment with the ACDM hypothesis, since the limiting mass be- 
low which all CDM-halos ar e dark ought to be two orders of 
magnitude lower according to lLi et al. (2010). 

As a final note, the newly derived reduced mass of Hercules 
(see end of Sect. 12.8b affects neither the calculated likelihoods 
nor the conclusions reached here. 

Summarising Sect. \3\ the mass distribution of the predicted 
DM halos of observed satellites is consistent with the ACDM 
hypothesis with at most 4.5 per cent likelihood. Assuming the 
dSph satellites are in virial equilibrium, the observationally de- 
duced DM halo masses of the MW satellites show a significant 
overabundance of Mo.3kpc ~ 1O 7 M halos and a lack of less- 
massive values compared to the theoretically calculated distri- 
bution for luminous sub-halos, despite much effort to solve the 
common-mass-scale problem (Sect. [2j. 



4. The bulge mass versus satellite number relation 
(problem 3) 

According to a straight forward interpretation of the CCM, 
more massive DM host halos have a larger number of luminous 
satellites because the number of sub-halos above a low-mass 
threshold increases with host halo mass, given the host halo 
mass waxes by accreting sub-halos. The sub-halos are accreted 
mostly individually without a physical link to the processes oc- 
curring at the centre of the host halo. There indeed does not ap- 
pear to be an observed relation between the halo mass and the 
bulge mass, si nce pairs of galaxies with and wit hout a bulge 
(such as M31, iRubin & Fordlll970l and M101, iBosma et al 



198 ll respectively) but with the same rotation velocity can be 
found. It would be useful to return to models A-F (Sect. [2]i 
and to include the formation of the host galaxy in the mod- 
elling to quantify the degree of correlation between the bulge 
mass and number of luminous satellites actually expected in 
the CCM. When doing so, the same type of models will also 
have to account for the presence of bulge-less galaxies hav- 
ing the same DM-halo mass, as pointed out above. That is, 
it would not suffice to merely demonstrate that some sort of 
bulge-mass-satellite number correlation emerges in the CCM. 
The case Mb u i ge = must emerge naturally within the model, 
since two-thirds of all bright dis k galaxies have no bulge or 
only a small one dCombesll2009bl) . 

On the basis of extragalactic observational data, 



Karachentsev et al.l (120051) note, but do not quantify, the 
existence of a correlation between the bulge luminosity and 
the number of associated satellite galaxies such that galaxies 
without a bulge have no k n own dSph companions, such 



The existence of this correlation can be tested in the Local 
Group, where good estimates of the number of satellites within 
the nominal virial radii of the respective hosts and of the stel- 
lar bulge masses of the three major galaxies (MW, M31, and 
M33) exist. Only the satellites brighter than L v = 0.2 x 10 6 L 
(My < -8.44) are considered, given that the census of fainter 
satellites is incomplete for the MW (notably in the southern 
hemisphere), and also for M31 and M33 given their distances. 
By restricting this analysis to satellites with Ly > 0.2 x 10 6 L o , 
the result becomes robust against the discovery of additional 
satellites since these would typically be fainter. The result is 
displayed in Fig. [3} a linear correlation between the bulge mass 
and the number of early-type satellites is suggested. An error- 
weighted least squares linear fit to the data yields 



N dSph = (4.03 ± 0.04) x M bulge /(1O 1O M ). 



(17) 



In terms of the present-day stellar mass fraction, the dSph satel- 
lites of the MW add-up to at most a few times 10 7 M , so that 
they amount to about 0.15 per cent of the mass of the bulge. 
Given that Eq. [T7] is a linear fit to three data points only, it 
will be important to check the reality of this correlation by sur- 
veying disc galaxies in the Local Volume with different bulge 
masses for a deep and exhaustive sampling of satellite galaxies. 

Given the small number of observational data points under- 
lying Eq. [T7J one should not over-interpret this result, but it is 
legitimate to inquire how significant the empirical correlation 
between bulge mas s and the number of satelli tes is. In view of 
the observation by iKarachentsev et al.l d2005l) noted above, it 
may be indicative of a physical correlation. 

The significance of the Local Group bulge-satellite corre- 
lation is evaluated by performing a Monte Carlo experiment, 
the null hypothesis of which is that there is no correlation. This 
hypothesis would appear to be plausible in the CCM because 
the number of satellites depends mostly on the host DM halo 
mass, while the bulge is produced by baryonic processes taking 
place near the center of the host DM halo. Three pairs of M bu i g e 
and A^asph values are chosen randomly from uniform distribu- 
tions such that M bulge e [0, 4.6 x 10 10 M Q ] and N dSph e [0, 280. 
For each three-point data set, a linear regression yields a mea- 
sure of the degree of correlation. This is performed 10 6 times. 
The following incidences are counted: 1) the resulting linear 
relation passes the (M bu i ge , A^dSph) = (0, 0) poinj| and the slope 
of the linear relation has a relative uncertainty smaller than a 
given value; and the second test is 2) the slope of the linear 
relation has a relative uncertainty smaller than a given value. 
The relative uncertainty in the slope used here is based on the 
uncertainties in the data. Applying this relative uncertainty to 
Eq.[l7]leads to N dSph ~ (4 ± 1) x M bulge /(1O 1O M ). Taking the 



as M101. IKarachentsev et al.l (120051) also point out that the 



number of known dSph satellites increases with the tidal 
environment. 



7 The upper bounds of the intervals are the 3cr upper values of M bu | ec 
and Misph of M31. The scaling of the axes is, however, irrelevant for 
the results of the Monte Carlo experiments, because the aim is to test 
how likely a correlation results, given the null hypothesis. 

8 The precise condition here is as follows: Let there be three Monte 
Carlo pairs (Mb u i ee , MiSph);, i = L..3. A linear regression yields a slope 
and an axis intersection, both with uncertainties expressed as <x val- 
ues. If the axis intersection lies within 5cr of the (0, 0) point, then this 
particular set of bulge-satellite pairs is counted. Note that the test does 
not require the slope to be the same as the observed value. 
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upper and the lower lcr limit of the slope, this equation thereby 
passes the lower and the upper lcr values of the data (Fig.[30 

The Monte Carlo result is that case 1) occurs 44 000 times, 
while case 2) occurs 157 000 times. Thus, if the correlation ev- 
ident in Fig. [3] were unphysical, then observing it would have 
a likelihood of 0.044 and 0.157, respectively. Given the data 
on the Local Group, the above hypothesis that the bulge mass 
and number of satellites are not correlated can therefore be dis- 
carded with a confidence of 95.6 per cent and 84.3 per cent in 
case 1) and case 2), respectively. 

Summarising Sect. the null hypothesis that the bulge 
mass and the number of satellites are independent quantities 
is rejected, based on the Local Group data, with a confidence 
of more than 95.6 per cent. With the absence of a DM-mass- 
luminosity relation for the observed satellites (Sect. [2]), this 
suggests that our present understanding of how satellite dwarf 
galaxies form and evolve may need revision. In the formation 
modelling of satellite galaxies within the CCM it will therefore 
be necessary to include also the formation of the host galaxy, to 
quantify the correlation between bulge mass and the number of 
satellites within the CCM. It will also be essential to refine this 
correlation using deep observational extra galactic surveys. 

5. The disc of satellites (DoS) and invariant 
baryonic galaxies (problems 4 and 5) 

The DoS is now addressed in view of the new satellite galaxies, 
and in Sect. I5.5l the issue that the two major DM halos of the 
Local Group, which happen to be similar, are occupied by sim- 
ilar disk galaxies is addressed within the context of the CCM. 

An important constraint for understanding the origin and 
nature of the observed satellite galaxies comes from them be- 
ing significantly anisotropically distributed about the MW, and 
possibly also about Andromeda. The pr oblem of the MW sys - 
tem for the CCM was emphasised by Kroupa et al.l d2005l) . 
They pointed out that the observed satellite system of the MW 
was incompatible at the 99.5 per cent confidence level with the 
theoretical distribution expected if the satellites were DM sub- 
halos tracing an isotropic DM host halo. Until then, the pre- 
diction within the DM hypothesis was that the distribution of 
sub-halos ought to be nearly sph erical and tracing th e shape 
of the host DM halo. For example. lAubert et al.l (120041) show a 
MW-type DM halo to have an infall asymmetry of only about 
15 per cent. The sub-halos enter the host halo along filaments 
and then phase-mix and virialise within the growing host DM 
halo. Similar sub-halo distributions are obtained in CDM and 



WDM models dKnebe et alj|2008l) . 

The DoS is a prono unced feature of the MW satellite sys- 
tem ( Metz et al. 2009bl) . an d a similar structure w as reported 
for the Andromeda system (IKoch & Grebell 120061) for which, 



9 The uncertainty in the slope given by Eq.[T7]is a measure for how 
close the data lie to the straight line fitted to them, i.e. very close in 
the given case. However, the uncertainties on the data suggests that 
the observed case is rather improbable (although obviously not im- 
possible), even if the correlation between A^sph and Mb u i gc is real. The 
uncertainty on the slope stated in Eq. fTT] would therefore not be a good 
basis for the test performed here. 




bulge 

Fig. 3. The number of dSph and dE satellite galaxies more 
luminous than 0.2 x 10 6 L o is plotte d vers us the bulge mass 



of the host galaxy (M W: IZhaoHl996[ M31: lKentl[l989l M33: 



Gebhardt et al.l 120011) . Only satellites within a distance of 



270 kpc of the MW and M3 1 are used. The solid line (slope= 
4.03) is Eq.[17] The upper (slope= 5.03) and the lower (slope= 
3.03) dotted lines illustrate the relative uncertainty assumed in 
the Monte Carlo experiment (see Sect. [4]). 

however, the distance uncertainties are larger and the satellite 
population is richer and more complex including dSph, dE, and 
dlrr galaxies. In the case of the well-studied MW, the DoS is 
very pronounced for the classical (11 brightest) satellites in- 
cluding the LMC and SMC. But how are the new satellites, 
the ultra-faint ones, distributed? Much hope for the CCM rests 
on the new discoveries being able perhaps to alleviate the DoS 

pro blem. 

Watki ns et al. 1 d2009h and Belokurov et al I d2010h reported 
the discovery of two new MW satellite galaxies, Pisces I 
and II, respectively, enlarging the total MW satellite system to 
24 satellites. Pisces I and II were found in the southern part of 
the SDSS survey area, making them the two first non-classical 
satellite galaxies found in the Southern Galactic hemisphere. 
Furthermore, distances to a number of the already known satel- 
lite galaxies have been updated in recen t years, most notabl y 



the new distance estimate for Boo II by IWalsh et al. I d2008l) 



which changes the distance from 60 to 42 kpc. An updated list 
of all currently known satellites is provided in Table. [2] upon 
wh ich the following a naly sis is based. 



Metz et al.l d2007l) and lMetz et alJ d2009al) employed a so- 
phisticated fitting routine to find the DoS. Here, an intuitive 
plane-fitting algorithm and a new disc-test are introduced. The 
plane-fitting algorithm leads to perfect agreement with the re- 
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Notes to the table: Data for the MW satellite s used for fitting the DoS. Se g 1 and 2 (marked in italics) are included in this list for reference, but they have not been included in the fitting 
because they appear to be diffuse star clusters dNiederste-Ostholt et al]|20"09T) . The positions are given both in Heliocentric coordinates (right ascension 122000, declination £2000, and Heliocentric 
distance r hc i io for epoch J2000.0) and in Galactocentric coordinates assuming the Sun to have a distance of 8.5 kpc from the MW centre. / MW gives the Galactic longitude with 0° pointing from 
the Galactic centre to the Sun. £>mw is the latitude as seen from the Galactic centre and tmw me radial distance from the centre of the MW. The coordinates were obtained using data from the 
references listed in the column labelled Ref., where more than one source is given, the distances to the satellites were obtained by error-weighted averaging over the available measurements. 
The satellite's line-of-sight velocities with respect to the Galactic standard of rest (GSR) are calculated assuming the Sun to move into the direction / = 90°, b = 0° (in Heliocentric, Galactic 
coordinates) with a velocity of either 220 kms -1 (1 'q 2 s ° r ) or 250 kmr 1 (Vq S ° r ). The measurement-uncertainties for the radial velocities reported in the respective papers (referred to in column Ref.) 
are reproduced in the column labelled Av. Finally, L v gives the satellite luminosities in the photo metric V-band; ag ain uncertainty-weighted averages are quoted when more than one reference 
is given in column Ref. Data marked with t have measured proper motions, listed in tab le 1 of ' 



the err or i s estimated to be 20 per cen t of the distance . Refere nce s: 1: Aden et alii 2009al) . 2 



Belokurov et al 



20091). 6: IBelokurov et all fcOlOh 7: IColeman et ail d2007bh 8: Ide Jong et al.l d2008l). 9: Ide Jong et al.l d201 



Metz et al.l d2008l), *: As no distance uncertain ties for Pisces II are availa ble in the literature, 
d2006l). 3: IBelokurov et ail d2007t). 4: IBelokurov etail d2008l). 5: 
L 10: IPairOra et al.l (120061). 11 : iGeha et alj d2009h . 12: iGreco et a 



-DlOh. _. 

Ibata et al.ld2006h. 14: lKoch et al.l d2009t). 15 : lKollmeier et alN2009h. 1 6: lKuehn et al.l d2008l). 17:lMartin et al.l d2008al). 18:lMartin et alJd2008bh. 19:lMateol dl998h. 20:lMoretti et al I Hop 



Belokurov et al.l 



Mufioz et alj d2006l). 22: iMufioz et al.i 201dh. 23: iMusella et al l d2009h. 24-lokamoto et aTu2008l). 25 : [Sakamoto & Hasegawal d2006l). 2 6- lsand et al.l d200S ) 27: ISand et al.l d201cK. 28: Isiegej 



d2006t). 29: ISimon & Gehal d2007t). 30: lwalsh et alj ( l2007t) . 3 1 : IWalsh et al. ll2008l) . 32: IWatkins et al. I d2009h . 33: IWillman et alj d2005al) . 34: IWillman eTaiTd2005bl) . 35: IZucker etail d2006al) 



36: IZucker etal.ld2006bl) . 37: lvan denBerglJdl99 
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suits obtained by Metz et al., and the new test allows an assess- 
ment of how discy the satellite distribution is. 

5.1. Parameters of the DoS 

A simple and straightforward method is described to calculate 
the DoS parameters /mw, £>mw, D p , anc ' ^> which are > respec- 
tively, the direction of the DoS-normal vector in Galactic lon- 
gitude and latitude, the smallest distance of the DoS plane to 
the Galactic centre, and the root-mean-square height (half the 
thickness) of the DoS. 

The positions of satellites on the sky and their radial dis- 
tances (compiled for convenience in Table |2]i are transformed 
into a Galactocentric, cartesian coordinate system assuming the 
distance of the Sun to the centre of the MW to be 8.5 kpc. The 
z-coordinate points into the direction of the Galactic North Pole 
and the Sun lies in the MW disk plane. 

The 3D coordinates are projected into two dimensions, 
plotting z against a projection onto a plane defined by the 
Galactic longitude /mw- This resembles a view of the MW 
satellite system as seen from infinity and from within the MW 
disc plane. The view of the satellite system is rotated in steps 
of 1 °. For each step, a linear fit is made to the projected satellite 
distribution. The linear fit is determined using the least squares 
method, demanding the satellite-distances, as measured per- 
pendicularly to the fitted line, to become minimal. This line 
constitutes a plane seen edge-on in the current projection. The 
two free parameters of the fit are the closest distance from the 
MW centre, Dp, and the inclination /?mw of the normal vector 
to the z-axis (a polar plane has /?mw = 0°). The plane-normal- 
vector's longitude is /mw, given by the projection. The fits are 
performed for each angle /mw between 0° and 360°. After half 
of a rotation, the view is identical to the one of 180° before, 
mirrored along the z-axis. 

For each angle /mw, the root mean square (RMS) height, 
A, of the satellite distribution around the fitted line is deter- 
mined. The normal vector to the best-fit disc solution (the DoS) 
to the full 3-dimensional distribution of the MW satellites is 
then given by those /mw and /?mw that give the smallest RMS 
height A min . 

To account for the uncertainties in the distance of the satel- 
lites, the major source of error, the procedure is repeated 1000 
times. Each time, the radial position of each satellite is ran- 
domly chosen from a normal distribution centered on the satel- 
lite's radial distance. It has a standard deviation given by the 
distance uncertainties to the satellite galaxies. Once a realisa- 
tion with varied radial distances is set up, the coordinate trans- 
formation into the Galactic coordinate system is performed. 
The parameters of the best fits are determined for each realisa- 
tion. Their final values are determined by averaging the results 
of all realisations, the standard deviations of their values are 
adopted as the uncertainties in the fits. 

Fitting all 24 currently known satellite galaxies within a 
distance of 254 kpc from the MW, the minimum disc height is 
found to be Amin = 28.9 + 0.6 kpc. This is more than 14<x away 
from the maximum height of A max = 55.7 ± 1.3 kpc obtained 
at a 90° projection of the data. Thus, the DoS is highly signifi- 



cant. The position of the minimum height gives the best-fit disc, 
the DoS. The normal vector defining the DoS points to /mw = 
156°.4± 1°.8 and has an inclination of /?mw = -2°.2 + 0°.6, i.e. 
is nearly perfectly polar. Dp, the closest distance of the DoS 
from the MW centre, is 8.2 + 1.0 kpc «: A min . 

5.2. A novel disc test 

Another test to determine whether the satellite galaxies are dis- 
tributed in a disc can be performed by comparing the number of 
satellites near the plane to the number further away: Let A; n be 
the number of all satellites that have a perpendicular distance of 
less than 1 .5 times the minimal disc height A m ; n from the line- 
fit. Accordingly, A ou t represents all satellites further away. Both 
A^in and A out are determined for each rotation angle, measuring 
distances from the line (i.e. plane viewed edge-on in the given 
projection) that fits the distribution in the given projection best. 
This is illustrated in Fig. [4] It shows an edge-on view of the 
best-fit plane, along with a view rotated by 90°. Both views see 
the disc of the MW edge-on. 

Figure [5] shows the ratio of galaxies found within the DoS 
to those outside (solid black line), % = Ni n /N ou{ . The situation 
is shown for the unvaried distances. If the MW satellites were 
distributed in a disc, % would approach a maximum when look- 
ing edge-on, while it will rapidly decrease once the projection 
is rotated out of the disc plane. It is a good test to discriminate 
a disc-like distribution from a spheroidal one. The latter would 
not lead to much variation in the ratio. 

It can be seen that % approaches a maximum close to the 
best-fit /mw- At the maximum, only two of the 24 satellite 
galaxies are found outside of the required distance from the 
disc. The maximum % is thus 1 1.0, situated only a few degrees 
away from the /mw that gives the smallest height. This has to be 
compared to the broad minimum of fi * 1 . The disc-signature 
is obvious, proving the existence of a DoS that incorporates the 
new satellites found in the SDSS. 

5.3. Classical versus new satellites: is there a DoS in 
both cases? 

In addition to the above analysis of all 24 known MW satel- 
lites, the analysis is also carried out separately for two distinct 
subsamples: the 1 1 classical, most-luminous satellite galaxies 
and the 13 new satellites discovered mostly in the SDSS. Each 
of them uses an own minimal height, given by the subsam- 
ple distribution, in determining %. If all satellite galaxies fol- 
low the same distribution, given by the DoS, a separate fit- 
ting should lead to similar parameters. If, on the other hand, 
the new (mostly ultra-faint) satellites follow a different dis- 
tribution, then this would challenge the existence of a DoS. 
It is worth emphasising that while the brightest satellites in a 
ACDM model of a M W-type halo may excep tionally form a 
weak disc-like structure (Libeskind et al \200% . none of the ex- 



isting CCM-based theoretical satellite distributions predict the 
whole luminous satellite population to be disc-like. 

Furthermore, comparing the results for the classical 11 
satellites with the ones obtained by the more sophisticated fit- 
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Fig. 4. Parameters of the MW DoS: the 3-D distribution of the 
MW satellite galaxies. The 1 1 classical satellites are shown as 
large (yellow) circles, the 13 new satellites are represented by 
the smaller (green) dots, whereby Pisces I and II are the two 
southern dots. The two open squares near the MW are Seg 1 
and 2; they are not included in the fit because they appear to 
be diffuse star clusters nearby the MW, but they do lie well in 
the DoS. The obscuration-region of +10° from the MW disc is 
given by the horizontal gray areas. In the centre, the MW disc 
orientation is shown by a short horizontal line, on which the po- 
sition of the Sun is given as a blue dot. The near-vertical solid 
line shows the best fit (seen edge-on) to the satellite distribu- 
tion at the given projection, the dashed lines define the region 
±1.5 x A min , A min being the RMS-height of the thinnest DoS 
(A m ; n = 28.9 kpc in both panels). Upper panel: an edge-on 
view of the DoS. Only three of the 24 satellites are outside of 
the dashed lines, giving N m = 21, iV out = 3 and thus a ratio 
of ft = Nfa/Nout = 7.0. Note the absence of satellites in large 
regions of the SDSS survey volume (upper left and right re- 



gions of the upper panel, see also fig. 1 in iMetz et al.l l2009a 



for the SDSS survey regions). Lower panel: a view rotated by 
90°, the DoS is seen face-on. Now, only 13 satellites are close 
to the best-fit line, 11 are outside, resulting in % — 1.2. Note 
that by symmetry the Southern Galactic hemisphere ought to 
contain about the same number of satellites as the Northern 
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Imw'vcx degrees 

Fig. 5. Testing for the existence of the DoS. The behaviour of 
% for each view of the MW, given by the Galactic longitude of 
the normal vector for each plane-fit. % = Ni n /N out is the ratio 
of the number of satellites within 1 .5 x A m i n (A m i n = 28.9 kpc), 
Ni n , to those further away from the best-fit line, N out , calculated 
for all 24 known satellites, as well as for the fits to the 1 1 classi- 
cal and the 13 new satellites separately (taking their respective 
RMS heights as the relevant Amin)- The disc-like distribution 
can be clearly seen as a strong peak close to Imw = 150°. Note 
that the position of the peaks are close to each other for both 
subsamples separately. This shows that the new satellite galax- 
ies independently define the same DoS as the classical satellite 
galaxies. 



ting technique used bv lMetz et al.l (12007b is a good test to check 
whether the present technique gives reliable results. 

The graphs for both subsamples are included in Fig. [5] the 
results for classical satellites are represented by the dashed yel- 
low, the new (SDSS) satellite galaxies by the dashed green line. 
Both are in good agreement not only with the combined sam- 
ple, but also with each other. They peak at their best-fit Imw, 
with each of them having an N out of only one galaxy at the 
peak. 

Applying the technique presented in Sect. 15 . 1 1 to calculate 
the DoS parameters, the new satellites have a best-fit disc with a 
normal vector pointing to /mw = 15 1°. 4+2°. 0, only five degrees 
away from the direction that was obtained by considering all 
known MW satellites. The inclination is £>mw = 9°.l ± 1°.0, 
again an almost perpendicular orientation of the DoS relative 
to the MW disc, being only 1 1 degrees away from the value 
determined before. The derived RMS height is A mm = 28.6 + 
0.5 kpc, essentially identical to the one given by all satellite 
galaxies. The minimum distance from the MW centre is Dp = 
18.3+ 1.3 kpc. 

The fitting to the 1 1 classical satellites leads to results that 
are in very good agreement, too. The best-fit position for the 1 1 
classical satellites is /mw = 157°.6 + l°.l andb MW = -12°. 0± 
0°.5, the height is found to be A = 18.3 ± 0.6 kpc, and the 
closest distance to the MW centre is Dp = 8.4 ± 0.6 kpc. This 
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is in excellent agreement with the results of iMetz et al.l (120071) . 
In that paper, the authors reported that Imw = 157°. 3, £>mw = 
-12°. 7, A m ; n = 18.5 kpc, and Dp = 8.3 kpc. This illustrates that 
the results are extremely accurate despite employing a more 
simple disc-finding technique. 

The agreement of the fit parameters for the two subsam- 
ples separately is impressive. Two populations of MW satel- 
lite galaxies (classical versus ultra-faint) with different discov- 
ery histories and methods define the same DoS. This shows 
that the new, faint satellites fall close to the known, classical, 
DoS (sDoS c i). Even without considering the classical satellite 
galaxies, the new satellites define a disc, DoS new , that has es- 
sentially the same parameters. This confirms the existence of a 
common DoS~DoS new ~DoS c i- 



5.4. The DoS - Discussion 

A pronounced DoS is therefore a physical feature of the MW 
system. But what is its origin? Is the existence of both the 
classical-satellite DoS c i and the new-satellite DoS new , such that 
DoSnew ~ DoS c i, consistent with the CCM? 

It has been suggested that the highly anisotropic spatial 
satellite distribution maps a highly prolate DM halo of the MW 
that would need to have its p rincipal axis ori ented nearly per- 
pendicularly to the MW disc (iHartwickKOOOh . However, there 
is still much uncertainty and disa greement as to the shap e and 



orientation of the MW DM halo: iFellhauer et all d2006l) used 
the bifurcation of the Sagittarius stream to constrain the shape 
of the MW DM halo to within about 60 kpc, finding it to be 
close to spherical. The measure ment of the shape of the DM 
halo of the MW within 60 kpc bv lLaw et al. I d2009l) . also based 
on the Sagittarius stream, suggests that the DM halo is triaxial, 
but with major and minor axes lying within the plane of the 
MW disc. The DM halo of the MW would therefore not trace a 
similar three-dimensional structure as the satellites, unless the 
major axis of the MW halo changes its orientation by about 
90 degrees beyon d 60 kpc and becomes e ssentially disc-like 
(i.e. highly oblate). Law & Maiewski ( 2010l) find a new slightly 
oblate solution to the MW DM halo shape at distances from 
20 to 60 kpc. In this solution, the minor axis points along the 
line Sun-MW-centre suggesting a similar orientation of this ex- 
tra potential as the DoS. The authors emphasise that this model 
is not strongly motivated within the current CDM paradigm, it 
merely serving as a "numerical crutch". Given this disagree- 
ment about the shape and orientation of the MW DM halo, a 
significant future observational and theoretical effort to clarify 
the situation is needed. 

An additional issue for the CCM is that the normal to the 
DoS is defined mostly by the outermost satellites, while the 
direction of the average orbital angular momentum vector is 
defined by the innermost satellites for which proper motions 
have been measured. Both, the normal and the average orbital 
angular momentum vector are nearly co-aligned implying a 
strong degree of phase-space correlation betwe en the satellites 
such that the DoS is rotating dMetz et al.ll2008l) . This rotational 
DoS is not expected if the satellites merely trace the MW DM 



halo, because they would have independent infall histories and 
would therefore not be correlated in phase space. 

This phase - space feature has been addressed by 



Libeskind et al 



d2009h . In a thorough analysis of struc- 
ture formation on MW-galaxy scales, they show that the 
MW constitutes an improbable but possible constellation of 
CDM-dominated satellites about a MW-type disk galaxy, the 
satellites having (of course) independent infall and accretion 
histories. 

They analyse an N-body sample of 30946 MW-mass DM 
host halos with mass in the range 2 x 10 11 M to 2 x 10 12 M 
for the properties of their substructure distribution. They first 
select from this sample only those halos that host a galaxy 
of similar luminosity as the MW (specifically, galaxies more 
luminous in the V-band than My = -20). From this remain- 
ing sample of 3201 (10 per cent) hosts, they select those that 
contain at least 11 luminous satellites, leaving 436 (1.4 per 
cent) host halos. In this sample of 436 systems, about 30 per 
cent have 6 luminous satellites with orbital angular momenta 
aligned to a degree similar to that of the MW system. Thus, 
only 0.4 per cent of all existing MW-mass CDM halos would 
host a MW-type galaxy with the right satellite spatial distri- 
bution. As the authors point out, this probability of 4 x 10" 3 
that the DM model accounts for the observed MW-type satellite 
system would be lower still if proper motion measurements of 
additional satellites a ffirm the orbital an gular momentum corre- 
lation highlighted bv lMetz et aD d2008l) . or if the satellites that 
may be discovered in the so uthern hemisphe re by the Stromlo 
Milky Way Satellite Survey dJerienll20l6Tr°l also lie within the 
DoS. All 13 new satellites define the same DoS as the 11 clas- 
sical ones, and furthermore, the latest additions in the southern 
Galactic hemisphere also lie in the DoS (Sect. 15.31 ), suggesting 
that the DM hypothesis is much less likely than 0.4 per cent 
to be able to account for the MW satellite system in MW-type 
DM halos. 



Li & H elmi (2008b and lD'Onghia & Lake! d2008l) propose 



an interesting alternative solution to the satellite phase-space 
correlation problem: they suggest that the correlation is caused 
by the infall of groups of DM-dominated dwarf galaxies. 
Unfortunately, this proposition is challenged by all known 
nearby groups of dwarf galaxies being spatially far too ex- 



tended to account for the thinness of the DoS (IMetz et al 



2009b). It may be thought that the groups that have fallen 
in correspond to compact dwarf groups that do not exist any 
longer because they have subsequently merged. But this is 
compromised by the observation that their putati ve merged 



counte rparts in the field do not s e em to exist dMetz et al 



2009bl) . Indeed. iKlimentowski etall d2010h model a MW-type 
system and deduce "... that such a disc is probably not an effect 
of a group infall unless it happened very recently" (their sec- 
tion 4.2.2). Furthermore, this notion would seem to imply that 
dwarf galaxy groups are full of dSph galaxies, while the pris- 
tine (before group infall) MW halo would have formed none, 
in conflict with the observed morphology-density relation (e.g. 
Okazaki & TaniguchfeOOOl) . 



http://www.mso.anu.edu.au/~jerjen/SMS_Survey.html 
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It needs to be emphasised that the DM-based models have 
so far not addressed the issue that the DoS lies nearly perpen- 
dicularly to the MW disc; DM-based models need to postu- 
late that this occurs, and it may indeed simply be chance. The 
combined probability that a DM-based model accounts for the 
observed MW-type system, which has the properties that the 
satellites have correlated angular momenta and form a DoS 
highly inclined to the baryonic disc of the host galaxy, can- 
not currently be assessed but is, by logical implication, smaller 
than 4 x KT 3 . 

But perhaps the MW is a very special system, an outlier 
within the DM-based framework? This possibility can be as- 
sessed by considering the nearest MW-similar DM halo. It 
hosts a similar disc galaxy, Andromeda, which has a similar 
satellite system to the MW but is however richer and more 
complex, and has a larger bulge m ass than the MW (Fig . |3). 
Andromeda may also have a D oS (IKoch & Grebelll2006l see 



also fig. 4 in lMetz et al J 12009 all I suggesting that these satel 



lite distributions may not be uncommon among MW-type DM 
halos. 

Thus, a Local Group consisting of two dominant DM halos 
of similar (MW-type) mass would have a likelihood of 0.4 per 
cent times 1.4 per cent, i.e. 5.6 x 10~ 5 , to appear with two MW- 
type disc galaxies, one of them having a pronounced rotating 
DoS with 1 1 or more luminous satellites, and the other having 
at least 1 1 luminous satellites. 

5.5. Invariant baryonic galaxies 

The Libeskin d et al. I d2009h analysis, described in Sect. 15.41 
also shows that about 10 per cent of MW-type DM halos 
would host a MW-luminous galaxy, the 90 per cent of oth- 
ers would presumably host galaxies with lower luminosities 
suggesting a large variation between DM halo and luminous 
galaxy properties. This however, appears to be a problem con- 
sidering the properties of observed disc galaxies. By using 
a p rincipal component analysis on hundreds of disc galax- 
ies, iDisnev et al.l (|2008) demonstrate that observed disc galax- 



ies are simple systems defined by one underlying parameter, 
rather than about six if the galaxies were immersed i n DM 
halo s. Citing additio nal similar results, Ivan den Berghl d2008l) 
and Gavazzil J2009h re ach the same con clusion, as well as 
Gentile et all d2009b and lMilgroml d2009al) . This is further sup- 
ported by an entirely independent study of star-forming galax- 
ies, whi ch again shows a remarkably small variation of be- 
haviour dPflamm-Altenburg & Kroupall2009bl) . The discovery 
that the ratio of DM mass to baryonic mass within the DM core 
radius is constant for galaxies (Sect. 16.4.11 below) is another 
statement of the same effect. 

The small amount of variation for disc galaxies thus ap- 
pears to be very difficult to reconcile with th e large variation 



inhere nt in the DM model, as quantified by the lLibeskind et al 



d2009h analysis: 90 per cent of MW-mass DM halos would have 
disc galaxies that differ substantially in luminosity from the 
MW in the CCM, and yet the closest neighbour, Andromeda, is 



similar to the MW. This is the invariant-baryonic-galaxy prob- 
lem. 

Summarising Sect. \5[ the CCM is highly significantly chal- 
lenged by the spatial distribution of MW satellite galaxies and 
by the similarity of rotationally supported galaxies. The phase- 
space correlation problem of the classical satellites is enhanced 
significantly after the inclusion of the new ultra-faint satellites, 
and the Local Group enhances the invariant baryonic galaxy 
problem. 

6. The origin of dSph and dE galaxies: The Fritz 
Zwicky Paradox, an alternative proposition and 
deeper implications 

What has been learned so far? The DM-mass-luminosity data 
of MW dSph satellite galaxies appear to be in conflict with 
the CCM results, and the mass function of DM masses of the 
dSph satellites is not in good agreement with the mass func- 
tion of luminous sub-halos calculated within the CCM. The 
correlation bulge-mass versus satellite-number is tentative (for 
having only three points) but will very likely pass the test of 
time because the error bars allow for a conclusive significance 
test. The two quantities appear to be physically related as indi- 
cated strongly by the Local Group data and also extragalacitc 
surveys, but clearly much more work needs to be done both 
observationally and theoretically to refine the implied corre- 
lation. The highly pronounced phase-space correlation of the 
MW satellites means that any formation mechanism must have 
involved correlated orbital angular momenta of the satellites. 

Given that the formation of a bulge involves highly dissi- 
pative processes, it emerges that a highly dissipative process 
seems to have formed the bulge and the correlated orbital an- 
gular momenta of the satellites. This leads to the real possibility 
that the origin of both the MW bulge and its satellite popula- 
tion is related to a galaxy-galaxy encounter. Indeed, it is well 
known and documented that galaxy encounters lead to the for- 
mation of bulges and tidal arms that can host the formation of 
tidal-dwarf galaxies (TDGs). These are then naturally corre- 
lated in phase space. Since the bulge and the satellites of the 
MW are about 1 1 Gyr old, we are led to the scenario that the 
proto-Galaxy may have had a major encounter or merger about 
1 1 Gyr ago during whi c h the bulge and the sa tellites formed 



( Pawlowski et al. 2010l) . Wetzstein et al. ( 20071) demonstrated 



1 1 Note that the rich satellite system of M31 may have a sub- 
population of satellites in a disc-like structure dMetz et al.l2009ah . 



in a series of numerical models that the number of TDGs in- 
creases indeed with the gas fraction of the pre-collision galaxy. 
This is relevant to galaxy encounters at a high redshift, where 
galaxy encounters are expected to have been frequent. 

Noteworthy is that a scenario for the origin of dSph satel- 
lite galaxies along the above lines had been suggested already 
before the DM hypothesis was widely accepted, namely that 
they m ay be ancient TDGs dLvnden-Bellll 19761 1 19831: iKunkell 
1979b . This proposition can naturally account for their corre- 
lated phase-space distribution in the form of a rotating disc- 
like distribution (Sect. [5}, and would lend a new viewpoint on 
the difficulty of understanding the properties of the MW dSph 
satellites as DM sub-halos doc umented above . Indeed, in a fa- 
mous conjecture, Fritz Zwicky dZwickvl 19561 on p. 369) states 
that new dwarf galaxies form when galaxies interact. As shown 
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here this leads to a contradiction with observational data when 
this conjecture is combined with his other famous conjecture 
according to which the mas ses of galaxies are dominated by 
Dark Matter dZwickvlll956l) . This contradiction is referred to 
as the Fritz Zwicky Paradox. 

6.1. The evolution of TDGs 

A natural way to explain the satellite phase-space correlation as 
well as the bulge-satellite relation is thus to identify the dSph 
satellite galaxies of the MW with a population of ancient TDGs 
that probably formed during a gas-rich encounter between the 
early MW and another galaxy. But if they all formed at the 
same time, how can the different chemical properties and star- 
formation histories of the different dwarf galaxies then be ex- 
plained within this scenario? If the DM hypothesis is not viable 
for the MW satellite population, how can the high mass-to-light 
ratios of the satellites be explained? 

It is known that the satellite gal axies all have ancient popu- 
lations of an indistinguishable age dGrebe]||2008l) . perhaps be- 
ing created when the TDGs were born. Or, the ancient popu- 
lation comes from the precursor host galaxy. TDGs may also 
form with globular clusters as long as the star-formatio n rate 
surpasses a few M G /yr for 10 Mvr dWeidner et al.H2004l) . The 
chemo-dynamical modelling bv lRecchi et alj ( 120071) has shown 



that once a TDG (without DM) forms it is not natural for it 
to blow out the gas rapidly. Rather, the rotationally-supported 
small gas-rich discs of young TDGs begin to evolve through 
self -regulated star formation either until their gas is consumed 
or removed through ram-pressure stripping. Consequently, 
their internal evolution through star formation can be slow and 
individual, such that TDGs that formed during one encounter 
event can exhibit different chemical properties many Gyr af- 
ter their formation. Removal of the interstellar medium from 
the TDG through ram-pressure takes about half to a few orbital 
times, which is typically one to a few Gyr after formation. This 
time scale is consistent with the obs erved cessatio n of star for- 
mation in most MW dSph satellites dGrebellll999h . The TDGs 
that have remained at large distanc es from their hosts retain 
their gas and appear as dlrr galaxies (IHunter et al Eoooh . Once 
formed, TDGs cannot fall back onto their hosts and merge since 
dynamical friction is insignificant for them. A TDG may be 
dispersed (but not accreted) if it happens to be on a near radial 
orbit, which, however, is unlikely given the torques acting on 
the tidally expelled material from which the TDG forms during 
the encounter. 

If the dSph satellites are ancient TDGs then un- 
derstanding their internal kinematics remains a challenge 



of DM dBarnes&He 


rnauist 1992: Wetzstein et al 2007 


Bournaud et all 20071 


Gentile et all 20071 Milsrom 2007b. 



However, the inferred large M/L ratios of dSph satellites (and 
especially of the ultra-faints) may not be physical values but 
may be a misinterpretation of the stellar phase-space distribu- 
tion within the satellite. If this were to be the case then the ab- 
sence of a "DM-mass"-luminosity relation (Sect. [2} for dSph 
satellites would be naturally understood. 



The following gedanken-experiment illustrates that this 
could be the case. An unbound population of stars on simi- 
lar orbits, each slightly inclined relative to the other orbits, will 
reconfigure at apogalacticon and an observer would see a stel- 
lar phase-space density enhancement and would also observe 
a velocity dispersion. The M/L ratio calculated from the ob- 
served velocity dispersion would not be a true physica l M/L 
ratio. Models related to this idea have been studied by Kuhnl 
(119931) . Moreover, resonant orbital coupling c an periodically 
inflat e kinematically m easured M/L values dKuhn & M iller 



19891: iKuhn et all 1 19961) . Fully self-consistent Newtonian N- 



body models have demonstrated that unphysically high M/ L ra- 
tios arise indeed if TDGs are allowed to orbit a host galaxy suf- 
ficiently long such that the remaining stellar population within 
the ancient TDG ad opts a highly non-isotropic phase-space 



distribution function dKroupal 11997c iKlessen & Kroupal 11998 



Met z & Kroupal 120071) . These models suggest that it may be 
wrong to use an observed velocity dispersion to calculate a 
mass for a dSph satellite. Thus, tidal shaping of TDGs over 
a Hubble time can produce remnant objects that have internal 
highly-anisotropic stellar phase-space distributions that would 
be falsely interpreted by an obser ver as correspon ding to a high 
M/L ratio, as explicitly shown bv Kroup"al ( 1997 ). Intriguingly, 
these models r eproduce the gross phy sical parameters of dSph 
satellites well ( Metz & Kroupa 2007b . and thus constitute the 
simplest available stellar dynamical solutions of dSph satellites 
constructed without fine-tuning. 

It is indeed remarkable how model RSI -5 of Kroupal 

dl997l) . shown here as a snapshot (Fig. [6]), is an essen- 
tially pe rfect match to the dSph satellite Hercules (see 
2 in 



fig 



Coleman et al. I2007al) discovered 10 years later by 



Belokurov et alj (120071) . The half-light radius is 180 pc in the 
model and 168 pc for Hercules, RS 1-5 has a ve locity dispersion 
of about 2.8 kms~' (table 2 in iRroupall 1997b . while Hercules 
has a measured v elocity dispersion of 3.72 ± 0.91 kins"' 



(lAden et al.l l2009al) . and the inferred mass-to-light ratio that 
one would deduce from velocity dispersion measurements 
based on the assumption of equilibrium is about 200 in both 
cases. Both RSI -5 and Hercules have luminosities agreeing 
within one order of magnitude (the model being the brighter 
one), yet RSI -5 has no DM. 

The TDG models for dSph sa t ellites proposed by 



Lynden-Bell 


(1976 


bv Kroupal ( 


1997b 



based on observed properties of TDGs, thus lead to a popu- 
lation of ancient TDGs that are in reasonable agreement with 
the obser ved luminosities, dimen sions, and M/L ratios of dSph 



satellites (M etz & Kroupall2007l) . These model-dSph satellites 



require no fine tuning of parameters but only assume the forma- 
tion about 10 Gyr ago of about 10 7 M heavy TDGs composed 
purely of baryonic matter. This theoretical framework of satel- 
lite galaxies does not imply any relation between luminosity 
and (wrongly inferred) "dynamical mass", in agreement with 
the lack of this relation (Sect. 0. And it would naturally ex- 
plain why the mass function of luminous DM sub-halos cannot 
account for the observations (Sect. [3). Within Newtonian dy- 
namics, this dynamical modelling over many orbits around the 
MW DM halo has demonstrated that even low-mass satellites 
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Fig. 6. Model RSI -5 from iKroupal dl997h (on the kpc grid) 
is plotted over the surf ace brightness contours of Hercules by 
Coleman et al.1 (I2007al) (celestial coordinate grid). The dashed 
and dotted curve are, respectively, the past and future orbit of 
RS1-5. 



do not easily disrupt unless they are on virtually radial orbits 
dKroupalll997llMetz & Kroupall2007t) . 

Summarising Subsect. \6.1\ the physics of TDG formation 
and evolution is sufficiently well understood to conclude that 
1) once formed at a sufficient distance from the host, TDGs 
will take an extremely long time to dissolve, if at all; and 2) 
the TDGs formed will naturally lead to a population of ancient 
TDGs that resemble dSph satellites. A bulge-mass-number of 
satellite correlation and a DoS arise naturally in this scenario. 



6.2. On the substructure problem 

The MW dSph satellites can therefore be understood as an- 
cient TDGs that formed within a DM universe. But on the 
other hand, the extensive modelling within the CCM strictly 
implies, if DM is cold or warm (but not hot), that MW- 
luminous galaxies must be accompanied by hundreds (with 
a slight dependence on the cold or warm nature of DM) 
of s hining albeit faint satellites, which are not of tidal ori- 



gin dKnebe et al.ll2008t IMaccio et a"1l2010t iBusha et alJ 



Koposov e t al. 2009). For example, 



Tollerud et al.l d2008l) con 



20101 



jecture that "there should be between 300 and 600 satellites 
within D = 400 kpc of the Sun that are brighter than the 
faintest known dwarf galaxies and that there may be as many 
as 1000, depending on assumptions." Deep follow-up obser- 
vations of the low S/N ultra-low-lum inosity satellite candi- 
dates introduced bv lWalsh et al. (2009) show that these are not 
dSphs as a population. These results show that there is not a 
significant number of missing, ultra-low-luminosity satellites 
{M v > -2, D < 40 kpc) in the SDSS footprint, i.e. an area cov- 
ering half of the Northern hemisphere (Jerjen et al., in prep.). 
This may be a problem because of the ACDM prediction that 
there should be a dozen additional satellites (My < 0, D < 



40 kp c) in a q uarter celestial spher e (e.g. fig. 4 in lKoposov et al 
20091: see also lCooper eTai1l2010l) . 

If the dSph satellites are ancient TDGs stemming from 
an early gas-rich encounter involving the proto-MW and 
probably contributing a collision product to the MW bulge 
(see Sect. [U, th en this would mean that the MW would 
have a severe substructure problem as there would not 
be any satellites with DM halos less massive than about 
10 10 M Q with s tars, in con fl ict wi th DM predictions pro - 
vided by, e . g.. iKnebe et al.l d2008h. IPiemand et alJ d2008„ 
krsha et al.1 J201dh . IMaccio et al.l d2010h . and iKoposov et alJ 
20091) . '. 'erhaps a few dSph satellites are ancient TDGs, such 
as the classical or nine brightest satellites, and the remainder 
are the DM dominated sub-halos? This possibility is unlikely, 
because the new satellites span the same DoS (Sect. 15.31 ) and 
because they do not form a population with physical properties 
that differ distinctl y from those of the classical satellites (e.g. 
Strigari et al]|2008l) . 

Summarising Subsect. \6.2\ based purely on the existence 
of the satellite phase-space correlation and the formation and 
survival of TDGs in a hierarchical structure formation frame- 
work the Fritz Zwicky Paradox emerges and the validity of the 
DM hypothesis must be questioned, because the dSph satel- 
lites cannot be two types of object at the same time, namely 
DM-dominated sub-structures and ancient DM-free TDGs. 



6.3. Early-type galaxies 

But if TDGs account for the dSph satellites of the MW, would 
they then not also be an important population in other envi- 
ronments? The_prodjic^ion_of_TDGsjn CCM has been cal- 
culated by Okazaki & Taniguchil ( 2000l) . Intriguingly, they find 
that TDGs naturally match the observed number of dE galax- 
ies in v arious environments. The result of lOkazaki & Taniguchi 
d2000l) is rather striking, since they find that within the CCM 
framework only one to two long-lived (i.e., bright) TDGs need 
to be produced on average per gas-dissipational encounter to 
cater for the population of dwarf elliptical (dE) galaxies and for 
the density-morphology relation in the field, in galaxy groups 
and in clustera 12 ! 

Viewing dE galaxies as old TDGs would be consistent 
with them deviating from the mass-radius, M(r), relation of 
pressure-supported (early-type) stellar systems. The dE and 
dSph galaxies follow a r oc M 1 ^ sequence reminiscent of tidal- 
field-dominated formation. All other pressure-supported galac- 
tic systems (elliptical galaxies, bulges, and ultra-compact dwarf 
galaxies) with stellar mass M > 1 6 M follow instead the r ela- 
tion r oc m°- 60±0 01 (see fig. 2 in lDabringhausen et alfeoOSl see 



12 Note that lOkazaki & Taniguchil j200Ct) write: "Adopting the 
galaxy interaction scenario proposed by Silk & Norman, we find that 
if only a few dwarf galaxies are formed in each galaxy collision, we 
are able to explain the observed morphology-density relations for both 
dwarf and giant galaxies in the field, groups of galaxies, and clusters 
of galaxies." They also state "The formation rate of TDGs is estimated 
to be ~ 1 - 2 in each galaxy interaction." and proceed to compare this 
number with the actually observed number of TDGs born in galaxy 
encou nters. This statement is at odds with the quotation in lBournaua 

J2oTch . 
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also f ig. 7 in lForbes et al J2008I and fig. 1 1 in lGraham & Worlevl 



20081), which may re sult from opacity-limited monolithic col- 
lapse (IMurravl 120091) . Viewing dE galaxies as TDGs would 



also be consistent with the observation that they have es 
senti ally stellar mass-to-light ratios similar to globular clus- 



ters dBender et all 1 9921 iGeha et alJ2003l: iDabringhausen et al 
l2008HForbes et al.ll2008l) . If dE (baryonic mass > 10 8 M@) and 
dSph (baryonic mass < 10 8 M ) galaxies are old TDGs, why do 
they appear as different objects? That the dE and dSph galaxies 
differ in terms of their baryonic-matter density may be a re- 
sult of the finding that below 10 8 M Q spheroidal objects on the 
r cc M 1 ^ relation ca nnot hold their warm gas and consequ ently 
they must expand ( Pflamm- Altenburg & Kroupal 2009a ). be- 
coming more susceptible to tides from their host. 

dE galaxies are pressure-supported stellar systems, whil e 
young TDGs are rotationally supported (Bo urnaud et al .120081) . 
With a mass of less than typically 10 9 M Q , the velocity disper- 
sion of their stellar populations becomes comparable to their 
rotational velocity (of the order of 30 kms -1 ). That a size- 
able fraction of dE galaxies show rotation, some even with spi- 
ral structure djerien et al. 2000; Barazza et al. 2002; Geha et al 
2003t iFerrarese et al.l 12006c IChilingarianl l2009t : iBeaslev et al 



2009), is thus also consistent with their origin as TDGs. For an 



excelle nt review on dE galaxies the reader is referred to lLisker 
d2009l) . 



One is thus led to the following seemingly logical impasse, 
i.e. to the Fritz Zwicky Paradox. In the CCM, TDGs are formed 
and their number and distribution is calculated to match the 
number and distribution of observed dE galaxies in the different 
environments. Within the CCM, the observed luminous dwarf 
sub-structures are thus naturally accounted for by TDGs. But 
the dE galaxies cannot be both, DM sub-halos and TDGs at the 
same time. 



Summarising Subsect. 16.31 the physical processes at play 
during structure formation in the CCM imply that dE galaxies 
ought to be identified as ancient TDGs. Thus, there would be 
no room for shining DM substructures. 

6.4. Deeper implications: gravitational dynamics 



In Sects. Kx2l and l6.3l it has been shown that the DM hypothesis 
leads to the Fritz Zwicky Pradox when accounting for the num- 
ber of satellite and dE galaxies because the formation of TDGs 
is an intrinsic outcome of structure formation. In Sects. |2]to[5]it 
has also been shown that the CCM seems to have a rather ma- 
jor problem accounting for the observed Galactic satellites and 
their internal properties. This situation suggests that alternative 
ideas should be considered to help us understand the origin of 
these problems, and indeed repeat the steps that had led to a 
full-fledged DM framework of structure formation but with a 
different outlook. Since structure formation in the DM frame- 
work relies on Newtonian gravitation in the weak-field limit, 
one is naturally led to relax insistence on Newtonian dynamics 
in the weak-field limit and to consider modified gravitation the- 
ories, which remain compatible with General Relativity in the 
strong field regime and with the observed large-scale structure. 
We note that adopting non-Newtonian dynamics in the weak- 



field limit would not necessarily rule out the existence of DM: 
on the scale of galaxy clusters DM might still be needed, but 



instea d of being warm or cold, it would be hot (I Angus et al 
120091) . 



6.4.1. Non-Newtonian weak-field gravity 

Alternatives to Newtonian dynamics in the weak-field 
limit have been studied in great detail. The increasingly 
popular modified-Newtonian-dynamics (MOND) approach 
rests on a modification of the Newtonian acceleration 
in the weak-field limit, i.e. when the N ewtonian accel 
eration a is smaller t h an a t hreshold an dMilgroml 1 1983 



Bekenstein & Milgroml 11984 ISanders & McGaughl 1 2002 : 



Bekensteinl2004 : Famaev & BinnevM2005 : Famaev et al.l 2007 



Sandersl l2007i 120081; iMcGaughl l2008t iNipoti et al.l 12008 
Tiret & Combes! 120081 iBruneton et al.l I2OO9I) . A modified 



gravity (MOG) adding a Yukawa-like force in the weak-field 
limit has also been un der investigation (IMoffat & Tothll2009at 
Moffat & Tofhl l2009bl and references therein). In addition, 



an extension of the General Theory of Relativity to a class 
of alternative theories of gravity without DM and based on 
generic functions f(R) of the Ricci scalar curvature R have 
been developed and successf ully applied to the prob lem of 



galactic rotation curves (e.g. ICapozziello et al. 120091) . For a 



brief review of MOND and MOG and Milgrom's proposition 
on the possible physical origin for the existence of ao, the 
reader is directed to the Appendix. 

Both the MOND and MOG approaches have been 
applied to the satellite galax y problem with apprecia 



ble success dMilgroml 1 19951 jBrada & Mil groml 



Angusl 120081; IMoffat & Tofhl 120081 iHernandez etal 



2000; 



2010; 



McGaugh & Wolf I l2010h . It has already been conclusively 



demonstrated that spiral galaxy rotation curves are well recov- 
ered in MOND purely by the baryon distribution without any 
parameter adjustments ( Sanders & McGau gq 20021 : iMcGaugh 



1011 

2004 l2005al ISanders & Noordermeerl 120071). and MOG is 



repor ted to also do well on this account dBrownstein & Moffat 



2006). In contrast, the DM approach can only poorly repro- 



duce the vast variety of rotation curves, a nd cannot explain 
the amazing regulari t ies found in them ( McGaugh 2004; 
McGaugh et ail 120071: iGentile et al] 120091; iMilgronj | 2009ab . 



Notably, the realisation ( Gentile et all 20091 Milgrom 



2009a) 



that the ratio of DM mass to baryonic mass within the DM 
core radius is constant despite the large variation in the 
DM-to-baryonic-matter ratio globally within galaxies cannot 
be understood within the DM hypothesis. A constant ratio 
within that radius implies that the distribution of baryonic 
matter is indistinguishab le from that o f the supposedly present 
DM (as already found by Bosmafl98~l ). This implies a hitherto 
not predicted near-exact coupling between DM and baryonic 
matter that does not arise naturally in the CCM, while outside 
that radius the effects of DM should become noticeable 
McGaughl d2010h . The only way to physically couple DM and 



baryons with each other to this degree would be by postulating 
the existence of an unknown dark force that acts only between 
DM particles and baryons. The modified DM cosmology 
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would then comprise inflation, dark matter, a dark force, and 
dark energy. 

In MOND models, this behaviour of gravity comes nat- 
urally. That the rotation curves would be purely defined by 
the baryonic matter distribution in non-DM models indeed 
would naturally ex plain the later finding based on a large sam - 

(12008 ' 



pie of galaxies bv | Disnev et al 



Gentile et al. (2009), 



and Milgroml ( 2009al) that disc galaxies appear to be governed 



by a single parameter. Furthermore, the high galaxy-cluster- 
galaxy -cluster velocities required to obtain the features of the 
Bullet cluster have been shown to be extremely unlikely in the 
CCM (Sect. JT), but these velocities are found to naturally occur 
in MOND (lAngus & McGaughll2008l) . Last but not least, the 
time-delay problem of the CCM mentioned in Sect. Q] would 
disappear naturally. 

6.4.2. A consistency check 

If it were true that the physical Universe is non-Newtonian in 
the weak-field limit, then a simple test would provide a con- 
sistency check: high dynamical mass-to-light ratios, (M/L)d yn , 
(derived assuming Newtonian dynamics) would not be due 
to DM but due to the dynamics being non-Newtonian in the 
weak-field limit and/or be due to objects being unbound non- 
equilibrium systems (Sect. 16.1b . Thus, taking MOND to be 
a proxy for non-Newtonian dynamics in the weak-field limit 
(MOND is, arguably, the currently available simplest alterna- 
tive to Newtonian dynamics in the weak-field limit), all sys- 
tems with non-stellar (M/L)d yn values (as derived in Newtonian 
gravity) would have to have internal accelerations roughly be- 
low the MONDian valut0 a B = 3.9 pc/Myr 2 . That is, all 
pressure-supported (spheroidal) stellar systems that appear to 
be dominated dynamically by DM would need to have an in- 
ternal acceleration a < a a . Note that the emphasis here is on 
pressure-supported systems since rotationally supported sys- 
tems have been extensively and successfully studied in non- 
Newtonian gravitational theories and because dSph and dE 
galaxies are mostly pressure-supported objects. 
Figure [7] shows the acceleration, 



M 0.5T- 

a(r e ) = G -z = G T 



(18) 



that a star inside a pressure-supported system experiences at 
the effective radius, r e , of its host system with luminosity 
spanning 10 4 to 10 12 L Q . Here M = 0.5 T L v is the stellar 
mass within r e and Ly is the absolute V-band luminosity in 
solar units. The stellar mass-to-light ratio in the V-band is 
T « 3 for collisionless systems (two-body relaxation time 
longer than a Hubble time), while T* 1.5 for collisional sys- 
tems, i.e. for systems that have evaporated a significant frac- 
tion o f their faint low-mass stars by means of energy equipar - 
tition ( Kruijssen & Lamers 2008 ; Kruijssen & Mieske 2009b . 
Values of (M/L)d yn as high as 10 can be expected for purely 
baryonic systems if these retain their stellar remnants and hot 



13 Note that this statement is approximately true for all non- 
Newtonian gravitational theories since they must account for the same 
non-Newtonian phenomena in the weak-field limit. 
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Fig. 7. Upper panel: The dynamical (M/L)d yn ratio (calcu- 
lated assuming Newtonian dynamics to be valid in the weak- 
field limit) in dependence of the luminos ity, Ly, for pressure 



suppo rted stellar systems following IDabringhausen et al 
(2008). Note that here dE (< 1O'°L ) and E (> 10 l() L o ) galax- 



ies are both plotted with the same symbol. Lower panel: The 
Newtonian acceleration (Eq. fT8l of a star located at the effec- 
tive radius within the host system in dependence of the host 
luminosity. The dashed line is «o- Note that M/Ld yn is high in 
pressure-supported stellar systems only when a < «o- In both 
panels: UCD=ultra compact dwarf galaxy. Comparing the up- 
per and lower panels shows that evidence ofDM(M/Ld yn > 10) 
appears only when a < ciq. 



gas. For example, the mass of an E galaxy may be comprised 
of only 30 per cent or less of stars, the rest consisting of 
stellar remnants and gas that cannot cool to form ne w stars 
(IParriott & Bregmanll2008t IDabringhausen et al.ll2009l) . mean- 
ing that T = 5 would be an underestimate in that case. Ultra- 
compact dwarf galaxies, UCDs (sometimes also understood as 
extremely massive star clusters), hav e high stellar M/L values 
perhaps due to a bot tom-heavy IMF (Mieske & Kroupal l2008) 
or a top-heavy IMF (Dabring hausen et al ■1 12009b . 

By comparing the two panels in Fig. [7] it is indeed evi- 
dent that only those systems with a < a„ show non-baryonic 
(M/L)d yn values. This is more clearly shown in Fig. [8] where 
the MOND prediction for the range of dynamical mass-to-light 
ratios measured by a Newtonist living in a MONDian universe 



20 



Kroupa et al.: Local Group tests of cosmology 



is plotted as a function of Newtonian acceleration. For this fig- 
ure, the MOND expectation for the mass-to-light ratio, which 
an observer who thinks to live in a Newtonian world would de- 
duce, was calculated as follows. Adopting a conservative value 
of the baryonic mass-to-light ratio Tb ar between 0.7 (for a glob- 
ular cluster with an old metal-poor population depleted in low- 
mass stars) and 5 (for an old metal-rich pop ulation), the predic- 
tion of MOND inside the effective radius is (IFamaev & Binnev 
20051; lAngus et alj|2009b 



(19) 



(M/L) dyn mond 

= 0.5 x T bar x (l + Vl + 4a /fl) . 

We note that, writing customarily x = g/a , where g is the 
actual full acceleration experienced by a ballistic particle (in 
MOND0 Eq. [p follows from the form of the transition 
MOND function lMilgromll 1 983b 



jx{x) = x/(l + x), 



(20) 



w hich is val i d up t o x « 10. The theoretical transition derived 
bv lMilgromld 19991) and mentioned in the Appendix would yield 
virtually the same result. 

The three classical dwarfs that lie outside the predicted 
MOND range for (M/L) dyn in Fig. [8] are UMa, Dra co, and 
UMi. UMa may have an anisotropic velocity dispersion (Angus 
20081) : Draco is known to be a long-standing problem for 
MOND, but the te chnique of interloper removal developed by 
Serra et al.1 (120091) could probably solve the problem, although 
this particular case remains open to debate; UMi is a typi- 
cal example of a possibly out-of-equilibrium system, as it is 
elongated with substructure and shows evidence of tidal tails 
(D. Martinez-Delgado, priv. communication). Ultra-faint dwarf 
spheroidals are expected to be increasingly affected by this 
kind of non-equilibrium dynamic s, as shown to be true even for 
Newtonian weak- field dynamics (Kro upal 19971 Sect. 16. U. and 
even more strongly so in MOND dMcGaugh & Wolf | |2010). 

Summarising Subsect. \K4\ well-developed non-Newtonian 
weak-field approaches exist and have been shown to account 
for galaxy properties much more succesfully than the CCM, 
which would need to be extended by a dark force to account 
for the observed strong coupling between DM and baryons. All 
known pressure-supported stellar systems ranging from ellipti- 
cal to dwarf satellite galaxies behave dynamically as expected 
in a MONDian universe. In DM cosmology, the association of 
highly non-stellar (M/L)d yn values with a < a a would be coin- 
cidental as it is not built into the theory. It is, however, natural 
in a MONDian universe for observers who interpret weak-field 
observations with Newtonian dynamics. 

7. Conclusions and perspectives 

We inhabit a Universe for which physicists seek mathemat- 
ical formulations. A successful formulation of gravitational 
physics, the General Theory of Relativity (GR), requires the ex- 
istence of non-baryonic dark matter (DM) in order to account 
for the observed rotation curves of galaxies and other dynami- 
cal effects in this theory, which has Newtonian dynamics as its 
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14 In the notation applied here, the MOND formula becomes a 
fi(x)g, where the Newtonian acceleration a is given by Eg .1181 



Fig. 8. The correlation between the acceleration a{r e ) and the 
dynamical mass-luminosity ratio (M/L)d yn derived assuming 
Newtonian dynamics is shown for the same objects as in Fig. [7] 
The shaded region indicates the range in (M/L)d yn as it fol- 
lows directly from MOND models (without any parameter ad- 
justments) using Eq. [19] The graph shows the consistency of 
the data in a MONDian universe for an observer who inter- 
prets observations with Newtonian dynamics. Encircled dwarf 
spheroidals outside this range (UMa, Dra, and UMi) may indi- 
cate non-equilibrium dynamics, either because the whole sys- 
tem is unbound, or because of unbound interloper stars among 
the member stars (see Sect. 16.4.21 ). That virtually all pressure- 
supported stellar systems fall in the shaded MOND region sug- 
gests a successful consistency check. That is, stellar dynamics 
is MONDian rather than Newtonian on galactic scales. 



weak-field limit. On the other hand, non-Newtonian weak-field 
gravitational theories have also been formulated to account for 
the "DM-effects" observed in galaxies. 

Finding a definitive test that distinguishes between these 
two different solutions to the problem of galactic dynamics 
and cosmological structure formation is difficult. Both DM and 
modified gravity are designed to solve similar problems, so the 
test must rely on subtle differences between the models and the 
observational data. Thus far, GR+DM+A+inflation (the CCM) 
accounts fo r the e mergence of structure on large scales, and 
Reves et al. (2010) were able to exclude certain versions of al- 
ternative gravitational theories that had a lready been known b y 
the respective community to be unstable ( Contaldi et alJE008 ). 
But, as shown here, the CCM appears to have insurmount- 
able problems on galaxy scales such that other alternative ap- 
proaches need to be studied. A speculative ansatz to perhaps 
solve the observed near-exact DM-baryon coupling in galax- 
ies within a DM-Newtonian cosmology would be to extend the 
CCM by postulating the existence of a dark force (DF) leading 
to a GR+DM+DF+A+inflation cosmology that should perhaps 
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be investigated in more detail in the future. The greatest dif- 
ferences between the two competing approaches (CCM versus 
non-Newtonian dynamics in the weak-field limit) are expected 
in the weak gravitational regime where the subtleties of non- 
Newtonian weak-field dynamics are most pronounced, which 
is why the constituents of the outer edges of galaxies allow the 
most stringent tests. 

This contribution has statistically assessed whether the ob- 
served properties of satellite galaxies in the Local Group, 
which are the result of structure formation in the weak-field 
limit, are consistent with the CCM. Given that a substantial 
number of independent research groups working in the tradi- 
tional CDM and WDM approaches have by now made firm 
statements about the dwarf satellite galaxies of the MW and 
Andromeda such that the missing satellite problem is deemd to 
be solved, the CCM can be further tested sensitively on these 
scales within the Local Group. 

Five new problems for the CCM on the scale of the Local 
Group and dwarf galaxies have been uncovered: (i) the ob- 
served absence of a mass-luminosity relation (Sect. [2] the DM- 
mass-luminosity problem); (ii) the mass function of luminous 
galactic satellites (Sect. [3] the mass function of luminous satel- 
lite problem); (iii) the observed relation between the bulge 
mass and the number of satellites (Sect. [4] the bulge-satellite 
correlation problem); (iv) the accordance with the Milky Way's 
disc-of-satellites of the recently detected ultra-faint dwarfs 
(Sect. [5] the phase-space correlation problem); and (v) the low 
probability that two neighbouring MW-type DM halos con- 
tain similar MW-type disk galaxies (Sect. 15.51 the invariant- 
baryonic- galaxy problem) . 

It is found that the CCM is consistent with the Local Group 
data with a combined probabilitvPl p <K 3xl(T 3 . The fiveprob- 
lems thus appear to rather strongly challenge the notion that 
the CCM successfully accounts for galactic structure in conflict 



with a vast volume of reported research (compare with iFanelli 



201(1 All these challenges constitute a strong motivation for 
numerous future observational and theoretical investigations. 
For instance, the disk of satellit es will have to be confi rmed 
by surveys such as Pan-Starrs ( Burgett & Kaiser 2009k an d 
the Stromlo Milky Way Satellite Survey (SMS) dJerierfeOlOh . 
Given the existence of the DoS and by symmetry, the Southern 
hemisphere ought to also contain about 16 satellites, such that 
the SMS survey is expected to discover about 8 new southern 
satellites (Fig.|4]i. It will also be essential to refine the correla- 
tion between bulge-mass and satellite-number with extragalac- 
tic surveys. On the theoretical side, more inclusive modelling 
is needed to address these challenges within the CCM while, at 
the same time, existing viable alternatives should be explored 
with more emphasis. 



15 Summarising the likelihoods, p, that the CCM accounts for the 
observed data in the Local Group are in the individual tests: (1) mass- 
luminosity data: p x < 0.3 per cent (Sec. [2}; (2) mass function of lumi- 
nous sub-halos: pi < 4.5 per cent (Sect.O; (3) bulge-satellite num- 
ber: p3 x 4.4 per cent (Sect.|4j; (4) a MW-type galaxy with at least 
11 satellites in a DoS: p^ = 0.4 per cent; (5) a M31-type galaxy with 
at least 11 satellites: p$ = 1.4 per cent (Sect, l5"l4l . Thus, the combined 
probability that the general CCM framework accounts for the Local 
Group is p <sc 3 x 10~ 3 . 



With this contribution, the following clues have emerged 
suggesting the need for a new scenario for the origin and na- 
ture of dSph satellite galaxies. The observed correlation be- 
tween bulge mass and number of satellites suggests that a link 
between these two quantities may exist. The phase-space cor- 
relation of the classical and ultra-faint satellite galaxies im- 
plies that angular momentum conservation played an impor- 
tant role in establishing the satellite distribution. Given that 
bulges form in dissipational encounters, during which angular- 
momentum conservation rearranges matter on Galactic scales 
to be in highly correlated phase-space structures (tidal arms), 
a natural path thus opens to understand the likely origin of 
satellite galaxies. Already in the 1970's a tidal origin for dwarf 
spheroidal galaxies was suggested, based on their arrangement 
around the Milky Way (Sect.[6]l. This solution does imply, how- 
ever, that the dSph galaxies are ancient TDGs and not DM sub- 
haloes. Furthermore, by logical implication, dE galaxies would 
also be TDGs (Sec. 16.3b . This would imply that the vast ma- 
jority of < 1O 1O M DM sub-halos are unable to make stars. 
This, however, would be in conflict with all the CCM computa- 
tions (the Fritz Zwicky Paradox) available to date to the extend 
that the CCM would have to be discarded in favour of a uni- 
verse without cold or warm DM. In this case, the non-Keplerian 
rotation curves of galaxies and other DM effects additionally 
suggest that we live in a non-Newtonian weak-field framework 
within which galaxies would be pure baryonic obiectj^l 

This scenario would naturally solve problems (iii) and (iv), 
while it would not imply a "dynamical mass"-luminosity re- 
lation if the dwarfs are out of equilibrium, so could possibly 
solve problem (i). For purely baryonic galaxies, problem (ii) 
would not exist anymore by definition. Problem (v) would also 
vanish naturally. What is more, while in the CCM the associa- 
tion of highly non-stellar (M/L)d yn values with a < a a would be 
coincidental because it is not built into the theory, it is natural 
in a non-Newtonian universe for weak-field observers who in- 
terpret observations with Newtonian dynamics. Noteworthy is 
that the same statement can be made for the Tu lly-Fisher scal- 
ing relation for rotationally-supported ga laxies dTullv & Fishei 
19771: lMcGaughll2005bl ICombesll2009aJ) as we ll as t he newly 



found s caling relation of iGentile et al.l ( 20091) and Milgroml 
d2009al) . The supposed mass-deficit seen in young rotating 
and gaseous TDGs (such as those of NGC 5291) constitutes 
independent empirical evidence towards this same statement. 
Young tidal dwarf galaxies (TDG), which should be devoid 
of collisionless DM, appear to nevertheless exhibit a mass- 
discrepancy in Newtonian dynamics. This is a significant prob- 
lem for the DM hypothesis , but it is natur ally explained by 



MOND (IGentile et al.l 12007; Mil gromll2007l) . Also, while the 



high Bullet-cluster velocity is hard to account for in the CCM, 
it is natural in MO ND (Sect [111641 and 16.4. j) . And, it has al- 
ready been noted by Sanders (119991) that the dynamical-mass - 



16 Given that Newton derived the gravitational 1 jr 2 law over a very 
limited physical range (Solar System), while with the Local Group 
gravitational physics is probed on a length scale nearly eight orders of 
magnitude larger and in a much weaker field regime, it need not be 
surprising that an adjusted gravitational law is needed. 
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Fig. 9. A new cosmological structure formation framework: 
the mangrove merger tree. In a modified-Newtonian weak-field 
framework, purely baryonic galaxies merge thereby spawning 
new dwarf galaxie s giving rise to the morphology-density rela- 
tion (adapted from lMetzll2008l) . 



baryon-mass discrepancy observed in galaxy clusters is nearly 
removed in MONDian dynamics. 

It would thus appear that within the non-Newtonian weak- 
held framework a much more complete, self-consistent, and 
indeed simpler understanding of the Galaxy's satellites as well 
as of major galaxies may be attained, than within the CCM. 

However, to affirm this statement, this alternative cosmo- 
logical scenario will have to be investigated in as much detail 
as is now available for the CCM in order to perform equivalent 
tests as presented here for the DM hypothesis and to ascertain 
which of the non-Newtonian weak-field dynamics theories (and 
which versions of the theories) can most successfully account 
for the physical world. Models of merging gas-rich disc galax- 
ies need to be computed in MOND, for example, to study how 
the formation of TDGs proceeds and how the number of satel- 
lites thus formed correlates with the bulge that forms as a result 
of the encounter. These populations of satellites associated with 
globular clusters that formed along with them would naturally 
appear in (more than one) closely r elated planes explaining the 
Lvnden-Bell & Lvnden-BeUd 19951) streams, because a gas-rich 
galaxy pair undergoes many close encounters in MOND, each 
spawning some TDGs and globular clusters, before perhaps fi- 
nally merging. 

Figure [9] schematically depicts the structure formation 
scenario in this non-Newtonian weak-field framework: while 
purely baryonic galaxies would merge, these events would 
spawn dwarf galaxies such that a density-morphology rela- 
tion woul d be established (more dE galaxies in denser envi- 
Okazaki & TaniguchilEoOOl). 



teraction times between galaxies are much longer, while the 
number of mergers is smaller than in a DM universe. Hence, 
the number of observed galaxy encounters would be given fore- 
most by the long time scale of merging and thus by more close 
galaxy-galaxy encounters per merging event rather than on a 
high number of mergers. 

This would imply that compact galaxy groups do not evolve 
statistically over more than a crossing time. In contrast, assum- 
ing DM-Newtonian dynamics to hold, the merging time scale 
would be about one crossing time because of dynamical fric- 
tion in the DM halos such that compact galaxy groups ought 
to undergo significant merging over a crossing time. The lack 
of significant evolution of compact groups, if verified observa- 
tionally, would appear not to be explainable if DM dominates 
galaxy dynamics. Analyses of w ell-studied compact g roups in- 
deed indicate this to be the case ( Presotto et alj|201ol) . 

Thus, many observational problems may be solved uncon- 
trived by adopting non-Newtonian weak-field dynamics, and 
perhaps this was, in the end, the most self evident expla- 
nation to the discovery of non-Keplerian rotation curves by 
Rubin & Fordl dl970lF . 
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Appendix: A brief review of MOND and MOG and 
Milgrom's proposition on the possible physical origin 
and value of ao 

Theoretical approaches tryin g to emb ed MO ND wi thin 



1 heoretical approaches tryin g to emb ed MO NU wi i 

Lorentz-covariant f ramework l|Bekenstein| 20041: SandersJ 



2005 



Zlosni k et al.l 120071: IZhad 120081: iB runeton & Esposito-Faresd 200S 
Blanche! & Le Tied 120091 : lEsposito-Faresd 120091 : ISkordiT I2OOS 
Milgromll2009bF are currently under intense scrutiny, and a quasi 



linear formulation of MOND has b een discovered only recently 
dMilgrorj201(il : IZhao & Famaey|201oh . which appears to allow easier 
access to N-body calculations. 

However, none of these theories is (y et) fully satisfactor 
from a fundamental point of view (see e.g. Contaldi et al.l 1200 

2oTch and 



IBruneton & Esposito-Fares e 2008; Reyes et all 



moreover 



none of them explains (yet) why the acceleration threshold, ao, which 
is the single parameter of MOND (adjusted by fitting to one single 
system), is about c VA/3 (where A is the cosmological constant and 
c the speed of light), or that ao » cHg/2n, where Ho is the current 
Hubble constant. They also require a transition function, fi(x) (e.g. 
Eq. I20t . from the Newtonian to the modified regime, a function not 
(yet) rooted in the theory. 

A possible explanation of the coincidence a ~ c VA/3 and 
a theo retically-based transition function are suggested by iMilgroml 
dl999h . In Minkowski (flat) space-time, an accelerated observer sees 
the vacuum as a thermal bath with a t emperature proportional to the 
observer's acceleration dUnrurJll975h . This means that the inertial 
force in Newton's second law can be defined to be proportional to 
the Unruh temperature. On the other hand, an accelerated observer in 
a de Sitter universe (curved with a p ositive cosmol ogical constant A) 
sees a non-linear combination of the Unruhl i 19751) vacuum radiation 



and of the iGibbon s & Hawking] dl977 ) radiati on due to the co smo 



logical horizon in the presence of a positive A. lMilgromlJl999h then 
defines inertia as a force driving such an observer back to equilibrium 
as regards the vacuum radiation (i.e. experiencing only the Gibbons- 
Hawking radiation seen by a non-accelerated observer). Observers ex- 
periencing a very small acceleration would thus see an Unruh radiation 
with a low temperature close to the Gibbons-Hawking one, meaning 
that the inertial resistance defined by the difference between the two 
radiation temperatures would be smaller than in Newtonian dynam- 
ics, and thus the corresponding accelera tion would be la rger. This is 
given precisely by the MOND formula o f lMilgroml ( ll983l) with a well- 
defined transition-function fi(x), and a = c(A/3) l/2 . Unfortunately, 
no covariant version (if at all possible) of this approach has been de- 
veloped yet. 

The theoretical basis of the MOG approach relies on chosen val- 
ues of integration constants in solving the equations of the theory. This 
approach seems to work well from an observational point of view, 
but it's fundamental basis needs further research, as is also the case 
for MOND. It is noteworthy th at a formulati on of MOG in terms of 
scalar, vector, and te nsor fields (lMoffaj|2006T) may possibly hint at a 



convergence with the lBekenstein d2004l) tensor-vector-scalar theory of 
gravity. 



